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Nomenclature 


ac 

aerodynamic center 

AR 

aspect ratio 

b 

wing span 

c 

chord 

a 

section lift coefficient 

C|a 

section lift slope 

CL 

aircraft lift coefficient 

Q 

section drag coefficient 

Cdo 

section zero lift drag coefficient 

Cb 

airrcraft drag coefficient 

Cdo 

aircraft zero lift drag coefficient 

CG. 

center of gravity 

QjOCV 

lift curve slope of vertical stabilizer 

Clow 

lift coefficient of wing 

CLcruise 

lift coefficient at cruise of aircraft 

CL5e 

change in lift due to elevator 
deflection 

Qna 

moment coefficient about C.G. 

CmO 

moment coefficient about C.G. at zero 
lift 

C n p 

yaw stability coefficient 

c root 

root chord 

e 

span efficiency factor 

E 

endurance 

Em 

max lift to drag ratio 

i 

current 

lx 

moment of inertia about x axis 

J 

propellor advance ratio 

Lia 

roll moment due to aileron deflection 

lv 

distance from C.G. to ac of vertical 
stabilizer 

mAh 

milliamp hours 

n 

load factor 

N.P. 

neutral point 

Q 

dynamic pressure 

R 

range 

RPM 

revolutions per minute 

S 

planform area 

S.M. 

static margin 

Sv 

planform area of vertical stabilizer 

t 

time 



X) 

take-off 

V 

volts 

V cruise 

cruise velocity 

V stall 

stall velocity 

Vto 

take-off velocity 

V v 

vertical tail volume ratio 

v„ 

freestream velocity 

a 

angle of attack 

Ocruise 

angle of attack at cruise 

a stall 

stall angle of attack 

r 

dihedral angle 

Secruise 

elevator deflection at cruise 

n 

propellor efficiency 

TW 

dynamic pressure ratios 


rolling angular acceleration 

<P 

bank angle 

e 

bank angle 

^max 

maximum allowable stress 

X 

control surface area ratio 

0) 

roll rate 



1.0 Executive Summary 


The Exodus Prime Mover is an overnight package delivery 

aircraft designed to serve the Northern Hemisphere of Aeroworld. The 
preliminary design goals originated from the desire to produce a large 
profit. The two main driving forces throughout the design process were 
first to reduce the construction man-hours by simplifying the aircraft 
design, thereby decreasing the total production cost of the aircraft. The 
second influential factor affecting the design was minimizing the fuel 

cost during cruise. The lowest fuel consumption occurs at a cruise 

velocity of 30 ft/s. Overall, it was necessary to balance the 1 economic 

benefits with the performance characteristics in order create a 
profitable product that meets all specified requirements and objectives. 

The SPICA airfoil section and a rectangular planform were 
selected to reduce construction hours neccessary to produce the wing. 
Its flat bottom and lift characteristics provide a balance between 
aircraft performance and construction simplicity. The wing area of 9.62 
square feet ensured the necessary lift both during cruise and take-off. 

In addition, cruise conditions occur at maximum lift to drag ratio. 

The Astro 15 electric motor and the ZingerJ 11-5 propeller 
comprise the propulsion system of the Prime Mover. The propeller 
selection was based upon the take-off distance requirement of 60 feet; 

the ZingerJ 11-5 provided the highest efficiency while still meeting this 
requirement. Twelve batteries of 1.2 volts and 1000 mah each were 
selected to power the system. The battery pack provides the voltage 
needed for take-off and the capacity required for the flight time of the 
aircraft. 

Directional and longitudinal control has been achieved through the 
use of a rudder and an elevator. A polyhedral concept has also been 
adopted for roll control. The polyhedral was chosen over the dihedral 
to decrease the amount of structure needed to withstand the bending 
moment at the root of the wing. 

The Prime Mover is capable of guaranteeing overnight delivery 
for the entire Northern Hemisphere due to the proposed fleet size of 42 
airplanes and the high range and endurance capabilities. The design 
objectives required the aircraft to meet a 8600 foot range minimum. 


The final design has displayed a cruise range of 24,000 feet, enabling 
the aircraft to complete its nightly schedule without the need to refuel. 
This reduces the operating costs of the aircraft. The maximum range 
and endurance of the fully loaded aircraft is 31,000 feet and 13.5 
minutes, respectively. The take-off distance at maximum take-off 
weight is 59 feet. 

The Prime Mover has a rectangular frontal area of 4.625 inches by 
4.375 inches and a fuselage length just under 5.0 feet to provide 800 
cubic inches of cargo space. The fuselage, wing, and empennage were 
designed to withstand a landing load factor of 4.0, a cruise load factor of 
2.5, and a catapult launch load factor of 2.0. 

The wing and the empennage will be removable in order to fit the 
disassembled aircraft within a 2ft x 2ft x 5ft box. Although this design 
increases the complexity of the structure, it enables the use of a 
modular construction technique. Each component of the aircraft may be 
built separately and assembled at a later time. This construction 
method will decrease the construction man-hours. 

As a result of the previously mentioned design characteristics. 
Exodus confidently presents the Prime Mover , an aircraft created to 
harmonize technical and economic considerations. The total production 
cost is estimated at $376,000. Based upon the production, operating, 
maintenance, and fuel costs Exodus recommends the price per cubic 
inch for intracontinental and overseas shipping be $8.74 and $11.01, 
respectively, in order to break even on the original investment. 







Prime Mover 


Specifications 


Aerodynamics 


Performance 


Wing Area 

9.62 ft 2 

Minimum Velocity 

10 ft/s 

Aspect Ratio 

9.62 

Maximum Velocity 

51 ft/s 

Chord 

12.0 in 

Stall Velocity 

22.3 ft/s 

Span 

9.62 ft 

Take-Off Velocity 

26.8 ft/s 

Taper Ratio 

1. 

Cruise Velocity 

30 ft/s 

Sweep 

0° 

R/C m ax 

5.4 ft/s 

Polyhedral 

17° 

Min Glide Angle 

3.33° 

Cd 0 (aircraft) 

0.018 

Take-Off Distance 

59 ft 

Airfoil Section 

SPICA 

Landing Distance 

44 ft 

Wing Mount Angle 

7° 

Max Range 

31,000 ft 

Stall Angle 

14° 

Max Endurance 

13.5 min 

CL max (aircraft) 

1.18 

Turning Radius 

60 ft 

L/D max (aircraft) 

17.8 



Structure 


Propulsion 


Empennage 


Engine 

Astro 15 

Hor. and Ver. Tail 

flat plate 

Propeller 

ZingerJl 1-5 

Airfoil Sections 


# of Batteries 

12 

Horiz. Tail Area 

208 in 2 

Total Voltage 

14.4 volts 

Vert. Tail Area 

77 in 2 

Battery Capacity 

1000 mahs 

Fuselage 




Length 

5.0 ft 



Height 

9.36 in 

Weights 


Volume 

1130 in 3 

Gross T.O. Weight 

6.7 lb 

Payload Volume 

800 in 3 

Empty Weight 

4.7 lb 

Frontal Area 

20.2 in 2 

Zero Fuel Weight 

3.45 lb 

Landing Gear 


C.G. position 


Type 

Tail 

x (from nose) 

30.5 in 


Dragger 

y (from bottom 

3.12 in 



of fuselage) 


Economics 




Fleet Life Cost 

$6.47 



/Volume Moved 




Overseas Price/in 3 

$11.01 



Intercont. Price 

$8.74 



/in 3 




Production Cost 

$376,000 





5 


Parameter 


Design Completion 


Prototype 

Completion 



DESIGN GOALS: 


V cruise 


Altitude cruise 


Turn radius 


Endurance 


Max Payload Volume 


Ranee-max payload 


Payload at Max R (w 


Ranee-min payload 


Desien life cycles 


Aircraft sales 


Tareet cost/in^ payload 


Tareet cost/oz payload 


30 ft/s 


25 ft 


60 ft 


13.2 min 


800 in 2 


24.000 ft 


0. oz 


26.200 ft 


6.7 lb 


600 


$376.000(actual cost 


$6.47(actual cost) 


N/A 


30 ft/s 


25 ft 


60 ft 


13.2 min 


800 in3 


24.000 ft 


0. oz 


26.200 ft 


6.8 lb 


600 


$376,000(actual cost 


$6.47(actual cost) 


N/A 





BASIC CONFIGURATION: 


Wing Area 


Weight (no 


Weight 


Wing loading (max Wgt) 


length 


span 


height 


width 


location of ref. axis origin 


nose/bottom 


9.62 ft 2 


4.7 lb 


6.7 lb 


0.7 lb/ft 2 


5.0 ft 


9.62 ft 


4.375 in 


4.625 in 


: fuselage 


6.8 


0.7 lb/ft- 


5. 


9.62 ft 


4.375 in 


4.625 in 


nose/bottom of fuselage 



WING: 


Aspect Ratio 


Span 


Area 


Root Chord 


Tip Chocrd 


Taper ratio 


-MAC 


Leading edge swee 


1/4 chord swee 


Dihedral 


Twist (washout 


Airfoil section 


Design Reynolds number 


t/c 


Incidence angle 


9.62 


9.62 ft 


9.62 ft 2 


12.0 in 


12.0 in 


1 . 


-.05 


0 . 


0 . 


17° 


0 . 


SPICA 


190,806 


11.7% 


9.62 


9.62 ft 


9.62 ft- 


12.0 in 


12.0 in 


1 . 


-.05 


0. 


17° 


0. 


SPICA 


190.806 


11.7% 


Hor. pos. of 1/4 MAC 

2.08% of chord 

2.08% of chord 

Ver. pos. of 1/4 MAC 

16.7% of chord 

16.7% of chord 












































































































c 


) - Wing 


C Lo - win ? 


- wing 


0.76 


0.0087 


0.1474 


0.0737/deg 



FUSELAGE 


Length 


Diameter - max 


Diameter - min 


Diameter - av 


Finess ratio 


Payload volume 


Total volume 


Planform area 


Frontal area 


CD 0 - fuselage 


C La - fuselage 


EMPENNAGE 


HORIZONTAL TAIL 


Area 


Span 


Aspect ratio 


Root Chord 


Tip Chord 


Taper ratio 


L.E. swee 


Incidence angle 


Hor. pos. of 1/4 MAC 


Ver. pos. of 1/4 MAC 


Airfoil section 


e - Oswald efficicnc 


Cd 0 * horizontal 


Clq ‘ horizontal 


-La * horizontal 


C Mm ac - horizontal 


VERTICAL TAIL 


Area 


Aspect ratio 


Root chord 


chord 


Taper ratio 


L.E. swee 


1/4 chord swcc 


5.0 ft 


N/A 


N/A 


N/A 


13.0 


800.0 in- 


1130.0 in3 


240.0 in 2 


20.2 in- 


0.004 


0.344 


208.0 in 2 


26.0 in 


3.25 


8.0 in 


8.0 in 


1 . 


0 . 


0 . 


+5° 


-2 ft from c 


+ 1.09 in from c 


late 


0.82 


0.0016 


0 . 


3.89/rad 


0.226/rad 


0 . 


0.76 


0.0087 


0.1474 


0.0737/dcg 



5.0 ft 


N/A 


N/A 


N/A 


13.0 


800.0 in 


1130.0 in3 


240.0 in 


20.2 in- 


0.004 


0.344 



208.0 in 


26.0 in 


3.25 


8.0 in 


8.0 in 


1. 


0. 


+5° 


-2 ft from c 


+1.09 in from c 


flat plate 


0.82 


.0016 



3.89/rad 


.226/rad 


























































































Hor. pos. of 1/4 MAC 


Ver. pos. of 1/4 MAC 


Airfoil section 


SUMMARY 

AERODYNAMICS: 


(airfoil) 


Lift curve slope (aircraft 


Cd 0 (aircraft) 


Alpha stall (aircraft 


Alpha zero lift (aircraft 


L/D max (aircraft 


Alpha L/D ma . (aircraft) 


-2.3 ft from c 


+5.76 in from c 


flat plate 


Avionics 


Pavload (max 


Motor 


Propeller 


Fuel 


Structure 


Win 


Fuselage/em 


Landin 


- max weight 


Number 


Placement 


Pavail max 


engine 


cruise 


Max. current draw 


Cruise current draw 


Propeller diameter 


Propeller pitch 


Number of blades 


. RPM 


Cruise Prop. RPM 


Max. Thrust 


Cruise Trust 


Number 


1.4 


1.18 


0.0737/de 


0.0179 


0.746 


14° 


-2 


17.80 


8.64° 


5.95 oz 


32 oz 


10.25 oz 


1.0 oz 


20.04 oz 


38.5 oz 


15.0 oz 


19.0 oz 


4.5 oz 


200 Watts 


15.4 Watts 


25 amps 


4.45 amps 


11.0 in 


5 in 


2 


18690 


4700 


1.5 lbs 


.38 lbs 


Panasonic 


12 


-2.3 ft from c 


+5.76 in from c 


flat plate 



0.0737/de 


0.0179 


.746 


14° 


-2 


17.80 


8.64° 



6.63 oz 


32 oz 


9.7 oz 


1.0 oz 


17.99 oz 


40.57 oz 


15.34 oz 


20.48 oz 


5.58 oz 



nose 


200 Watts 


15.4 Watts 


25 amps 


4.45 amps 


11 . 


5 in 


18690 


4700 


1.5 lbs 


.38 lbs 


Panasonic 


12 














































































Individual capacit' 


Individual voltage 


■H 


Pack voltage 


1000 mahs 


1.2 volts 


1000 mahs 


14.4 volts 


1000 mahs 


1.2 volts 


mahs 


14.4 volts 



STABILITY 

AND 

CONTROL: 



Neutral 


Static margin %MAC 


Hor. tail volume ratio 


Ver. tail volume ratio 


Elevator area 


Elevator max deflection 


Rudder Area 


Rudder max deflection 


Aileron Area 


Aileron max deflection 


42.4% chord from LE 


15.6% 


0.3 


0.01 


39 in 


- 20 ° 


49.5 in2 


+/- 20 


N/A 


N/A 


-0.672/rad 




3. 89/rad 




Endurance @ R max 


Endurance max - 


Range @ E^ 


Min Glide angle 


T/O distance 


T/O rotation angle 


Landing distance 


Catapult Range 


10.0 ft/s 


51.0 ft/s 


22.3 ft/s 


31,000 ft 


11.5 min 


13.5 min 


27,000 ft 


5.4 ft/s 


3.33° 


59.0 ft 


11.5° 


44.0 ft. 


not available 




SYSTEMS: 


Landin 


Main gear 


Main gear length 


Main gear tire size 


tail gear oostion 


tail gear length 


tail gear tire size 


Engine speed control 


Control surfaces 


tail dragger 


6.0 in 


5.0 in 


1.25 in 


57 in 


4.7 in 


1.1 in 


Tckin Speed Controller 


elevator/rudder 


42.4% chord from LE 


15.6% 


.01 


39 in 


+/- 20 ° 


49.5 in2 


+/- 20 ° 


N/A 


N/A 


-0. 672/rad 


.054/rad 


3. 89/rad 


.226/rad 



10.0 ft/s 


51.0 ft/s 


22.3 ft/s 


31,000 ft 


11.5 min 


13.5 min 


27,000 ft 


5.4 ft/s 


3.33° 


59.0 ft 


11.5° 


44.0 ft. 


not available 



tail dragger 


16.25 in 


5.0 in 


1.25 in 


57 in 


4.7 in 


1.1 in 


Tekin Speed Controller 


elevator/rudder 


















































































































TECH. DEMONSTRATOR: 


Payload volume 


Pavload weight 


Gross Take-Off Weight 


Weight 


Zero Fuel Weight 


Wing Area 


Hor. Tail Area 


Ver. Tail Area 


C.G. position 


1/4 MAC position 


Static margin %MAC 


Range max 


Endurance max 


V cruise 


Turn radius 


Airfram Struct. Weight 


Propulsion Syst. Weight 


Avionics Weight 


Landing Gear Weight 


Est. Catapult Range 


800.0 in 


32.0 oz 


6.7 lb 


4.7 lb 


3.45 lb 


9.62 ft 2 


208.0 in 


77.0 in 2 


30.5 in 


30.25 in 


17% 


26.8 ft/s 


31.000 ft 


13.5 min 


30.0 ft/s 


60.0 ft 


38.5 oz 


11.25 oz 


5.95 oz 


4.5 oz 


not available 


800.0 in3 


32.0 oz 


6.8 lb 


4.8 lb 


3.67 lb 


9.62 ft 2 


208.0 in 


77.0 in 


30.5 in 


20.25 in 


15.6% 


26.8 ft/s 


31.000 ft 


13.5 min 


30.0 ft/s 


60.0 ft 


40.57 oz 


10.7 oz 


6.63 oz 


5.58 oz 


not available 


ECONOMICS: 


Unit materials cost 


Unit propulsion system cost 


Unit control system cost 


Unit total cost 


Scaled unit total cost 


Unit production manhours 


Scaled production costs 


Total unit cost 


Single Flight 


Single flight o 


Single flight 


ross income 


. costs 


rofit 


#flights for break even 


$ 127 


$213 


$300 


$640 


$256,000 


120 


$ 120,000 


$376,000 


$6.47 


$2891 


$1638 


$1253 


12.600 


$164 


$213 


$300 


$ 677 


$270,800 


112 


$112,000 


$392,800 


$6.51 


$2891 


$1638 


$1253 


12.600 
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2.0 Mission 

2.1 Design Requirements and Objectives 

After carefully reading the request for proposal. Group Exodus 
strove towards developing a mission that would result in the maximum 
return on investment. In order to do this, the design requirements had 
to be assessed. Some of the primary requirements were as follows. The 
aircraft had to be able to take off and land under its own power within 
75 ft. . It had to be able to sustain a level 60 ft. radius turn. The 
aircraft had to travel at a velocity less than 30 ft. /sec. The 

disassembled aircraft had to fit within a 2’x2’x5’ box. The aircraft had 
to be able to reroute to the nearest airport and loiter for 1 minute in 
case of an emergency. Lastly, the aircraft could utilize no more than 4 
servos. These were thought to be the most important design 
requirements. 

Group Exodus was now prepared to analyze the Aeroworld market 
and develop the most profitable mission. In doing so, a set of design 
objectives was formulated. The complete list of these objectives is as 
follows: 

Structural: 

1. The aircraft design should be as simple as possible in order to 

minimize construction costs. 

2. The aircraft should maintain structural integrity throughout a 
life cycle of 600 flights. 

3. The aircraft must be able to survive a forced landing. 

4. The total payload volume will be between 576 and 768 cubic 

inches. This will allow for a minimum frontal area while 
maintaining a fuselage length of under five feet. 

5. Battery exchange should take less than 1 minute to reduce 

maintenance costs. 

6. The cargo must be easily accessible to decrease handling 
time. 



Performance: 

1. The average speed must be greater than 16 ft./sec so that 

the furthest destination can be reached in the allotted time. 

2. Aircraft must complete the furthest run twice within 12 

minutes to allow time for package handling. 

3. The maximum aircraft weight will be less than 7 pounds, the 

perceived maximum weight for previous aircraft of this class. 

4. The minimum range must be greater than 8600 ft. in order to 

accommodate the furthest projected flight. 


Economic: 

1. The total cost of the aircraft structure should not exceed 

$ 200 . 00 . 

2. Based on initial estimated production cost of $340,000 per 
aircraft, which is believed to comprise of 90% of the total 

life cost, the average price per cubic inch will be 
approximately $2.25. 

There were a few alterations to our original design objectives and 
these will be discussed in the following chapters. These objectives 
would now guide Group Exodus in formulating the complete mission 
statement. 

2.2 Mission Analysis 

Group Exodus, in arriving at the mission, realized the 

importance of the fact that we were designing for an inexperienced 
company in a new and uncharted market. Consequently, the fleet 
was designed for a restricted area which would result in excellent 
efficiency and dependability; as opposed to spreading the fleet thin 
throughout all of Aeroworld and decreasing the efficiency and 
dependability of the service. This dependable service would 
consequently attract more customers. Therefore, the Prime Mover 
fleet was designed to operate in the Northern Hemisphere only. 

There was another good reason to stay within the Northern 
Hemisphere. After further analyzing the economics of the mission. 



fuel costs were estimated to be approximately 73% of the total fleet 
life cost. Therefore, restraining the total mileage flown by the fleet 
would result in a less expensive service for the customer. As a 
result, a hub system was decided on. By centrally locating the hub 
city, the fuel cost for each plane could be kept to a minimum 
consequently reducing the total fleet life cost. City J was chosen as 
the hub city due to its central location and the large daily cargo load 
of this city. 

Another reason for chosing the Northern Hemisphere was the 

results calculated from the market analysis. After reviewing the 
daily cargo load departing from all of the Southern hemisphere cities 
(C, D, E, and O), there was no city that was delivering more than 1000 
cubic inches. City C ships 980 cubic inches, city D ships 690 cubic 
inches, city E ships 960 cubic inches and cty O ships only 280 cubic 
inches. After refering to Appendix I, it can be seen that the Northern 
Hemisphere city with the least amount of cargo, city H, ships almost 
590 cubic inches more than the Southern Hemisphere city that ships 
the most, city C. The Southern Hemisphere accounts for only 9% of 
the total departing cargo load of the world. Consequently, besides 

expending a lot of money in fuel costs, the prospect for profit in the 
Southern Hemisphere is quite weak and therefore Group Exodus 

recommends that only the Northern Hemisphere should be served at 
this point. Though our analysis does not include the Southern 
Hemisphere, the aircraft does have the capability to accommadate 

this area. Future servicing of this hemisphere may be considered 
using a different pricing system. 

There will be 42 planes in the Prime Mover fleet. The cargo 
area of each plane will be 800 cubic inches. This exceeded the 768 
cubic inch objective which was set to allow for single-stacked cargo. 
Single-stacked cargo reduces the frontal area and allows for better 
control over the cargo CG. We were able to increase our allowable 
cargo capacity by making the empennage removable, thus enabling 
the fuselage to be longer without exceeding the five foot storage 
restriction. This increase in volume will reduce the number of planes 
in the fleet and consequently greatly reduce the total fleet life cost. 
Each plane will make only two flights per day resulting in low fuel 



cost per plane as mentioned before. The plane will simply fly to city 
J, unload its cargo, reload, and return to its original destination. The 
planes that arrive at the hub city first can be unloaded and then the 
other planes can be unloaded as they arrive. Therefore, the loading 
and unloading process runs smoothly and is accomplished in less 
time. There is also the possibility of rotating the planes that fly a 
longer distance with those that fly a relatively short distance. This 
will even out the effects of fatigue on each of the aircraft. The 
simplicity of this mission will also result in basically no problems 
with the designing of routes and the scheduling of flights. This 
service will be able to accommodate 100% of the daily cargo load in 
the Northern Hemisphere or Aeroworld. The Prime Mover was also 
designed to be able to handle a full load of parcels at the maximum 
density. As a result, no customer will be turned away. Creating this 
peaceful state of mind in our customers will result in more customers 
which translates into more of a profit. 

The longest flight time for any route will be no more than 9 
minutes. This includes a flight to the hub city and back with a one 
minute loiter during each flight. This will leave 6 minutes for 
loading, unloading, and daily servicing. This is more than enough 
time in Aeroworld. As a result of this efficient time schedule, the 
service becomes as dependable as possible. The price of one aircraft 
will be $376,000, in order for the company to break even in half of 
the fleet life. The price breakdown is available in Section 12.1. 



3.0 Concept Selection Studies 


The design process was begun by developing a number of 
design concepts to meet the request for proposal. Acknowledging the 
strengths and weaknesses of each proposal allowed Exodus to 
determine the best possible final aircraft design. 

3.1 Concept #1 

Concept one is a conventional aircraft with a high mounted 
rectangular planform wing with dihedral. It has a conventional tail 
and a single forward mounted engine. This concept emphasizes 
simplicity, reducing construction hours and thus construction costs. 

Control would be achieved by use of the elevator and the 
rudder-dihedral combination. Therefore, only two servos would be 
required for the control surfaces, which decreases operational costs. 
The dihedral joint would be subject to a large moment due to it s 
location at the root of the wing. This would require a costly increase 
in structural support. 

Tail dragging landing gear is used to maintain landing stability. 
It allows for the option of attaching the rear landing gear to the 
rudder for increased ground control. 

3.2 Concept #2 

Concept two is a canard configuration. The main wing is high 
mounted with dihedral and the canard is mounted low and forward 
on the fuselage. This configuration of the lifting surfaces would 
minimize the inerference of the control surfaces on the main wing. 
Again the aircraft would be driven by one forward mounted engine. 
For this design, tricycle landing gear would implemented. 

A canard configuration would require the construction of two 
wings; hence, increasing the production hours and cost. The extra lift 
produced by the canard could be used to lift a larger payload. It 



becomes necessary to consider the power required to take off with 
such a large weight. 

3.3 Concept #3 

Concept three uses a low mounted wing with a tapered 
planform. The empennage is a T-tail, and the landing gear is a tail 
dragger. 

The tapered wing decreases the induced drag by increasing the 
wing efficiency and aspect ratio. However, the spanwise change in 
airfoil sections would increase the construction time and cost. 
Though the use of a T-tail would be beneficial in avoiding downwash, 
construction simplicity would suggest placing the horizontal stabilizer 
on the top of the fuselage Structurally, the empennage of a T-tail 
would be subjected to larger bending moments than a conventional 
tail. Control difficulties would arise in running the control wire from 
the fuselage to the top of the vertical stabilizer. 













Table 3.1 
Concept Evaluation 


Concepts Advantages Disadvantages 


Concept #1 

-Reduced 

Construction 

Costs 

-Landing Stability 
-Less 

Operational Cost 

-Increased Drag 
-Increased Wing 
Weight 

Concept #2 

-Extra Lifing 

Surfaces 
-Decreased 
Downwash 
Effects 

-Increased Weight 
-Increased 

Construction 

Costs 

Concept #3 

-Decreased 
Downwash 
-Decreased Drag 

-Increased 

Construction 

Costs 

-Empennage 

Structural&Control 

Difficulties 


3.4 Final Concept 

The Prime Mover took on a concept similar to that of the first 
concept. The wing is mounted above the fuselage in order to 
maintain a continuous cargo area. The high mounted wing is 
inherently stable and reduces the size of dihedral angle necessary to 
control the aircraft. A larger dihedral angle must be complemented 
with a larger wing area to produce the equivalent lift. This would 
result in a heavier wing. However, the actual area increase is fairly 
small. 

Exodus decided to use a polyhedral wing in place of a dihedral 
wing. The advantages of the dihedral are used without subjecting 




the wing juncture to the high bending moments located at the root of 
the wing. Lift of the Exodus Prime Mover is made comparable to that 
of a tandem design by using a large wing area. 

The wing planform is rectangular in order to decrease the 
construction hours and cost. For the same reason, the empennage 
will consist of flat plate stabilizers located on top of the aircraft. 

The Prime Mover will implement a single forward mounted 
engine and tail dragger landing gear. Exodus feels this design will 
fully meet the requirement and objectives of a profitable overnight 
delivery service. 


4.0 Aerodynamic Design Detail 


4.1 Airfoil Selection 

The process of selecting an airfoil began with the determination 
that the airfoil must demonstrate good characteristics at low 
Reynolds numbers in the estimated range of 150,000 to 200,000. 
Airfoils in this category were then judged on three design criteria: 
lift, drag, and geometry of the airfoil. An optimal airfoil design for 
our aircraft has a high Cimax at a large stall angle and gradual stall 
characteristics, minimum drag, flat bottom and thick trailing edge. 
The geometric characteristics were identified as important because 
they result in less construction time, and therefore reduce production 
costs. Four airfoils that best fulfilled the design criteria were 
selected from reference 7 for closer examination. These airfoils 

included the Wortmann FX63-137, Clark-Y, S3010, and the SPICA. 
The merits of each airfoil are listed in Table 2.1. 


Merits of Selected Airfoils 


Airfoils 

Merits 

Wortmann FX63-137 

High Cimax, High lift curve 

Clark-Y 

Thick trailing edge 

S3010 

Low Cdo 

SPICA 

Large a sta ll. Flat bottom 


The Wortmann FX63-137 has the highest Cimax of 1.6 at a stall 
angle of 12 degrees but the cambered airfoil with the thin trailing 
edge would be difficult to duplicate and especially difficult to 
monokote, thus greatly increasing production costs. The Clark-Y has 
a relatively flat bottom and the thickest trailing edge of the four 
airfoils examined. However, the airfoil's lift characteristics are poor 
for it has a Cimax of only 1.2 and stalls at an angle of 10 degrees. The 
S3010 airfoil has lower drag at zero lift than the other airfoils with 
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lift characteristics similar to the Clark- Y. Also, the bottom is not as 
flat as the Clark-Y and the trailing edge is thin. 

The SPICA airfoil was selected primarily for economic reasons. 
Its flat bottom and relatively thick trailing edge as shown in Figure 
4.1 make machining and handling easier, and therefore reduce 
production costs. In addition, the airfoil produces high lift 

coefficients with a Ci ma x of 1.45 at a Reynolds number of 151,000 

and stalls at an angle of 14 degrees as illustrated in Figure 4.2. The 
large stall angle is important so that the wing may be mounted at an 
incidence angle for increased lift at cruise with an adequate margin 
of safety. The airfoil exhibits gentle stall characteristics. However, 
the airfoil has high drag at low lift coefficients with a C<jo = 03. The 
SPICA airfoil was selected despite this fact due to the early 
determination that production costs were 90% of the total cost of an 

aircraft. It is important to note that greater consideration would 

have been given to the high section drag coefficients had the actual 
73% fuel costs been realized earlier. 

Figure 4.1 
SPICA Airfoil 



Thickness 11.70% 


Camber 4.75% 



4.2 Wing Design 

Once the airfoil was selected, the aerodynamic analysis focused 
on the design of the wing. The goal was to design a cost effective 
wing that provided sufficient lift during take-off, cruise, and landing 
with minimum drag. The design parameters for the wing included 
the planform area, aspect ratio, taper, twist, and sweep. 

Difficulties in constructing a tapered wing and a wing with 
twist would greatly increase construction time and costs as well as 
the possibility of wing warp, and therefore were not considered in 
the design process. Also, sweep was not considered since the aircraft 
was designed to fly at a low cruise speed of 30 ft/s and decreasing 
the effective velocity along the wing was undesirable. Therefore, a 
rectangular planform was selected for the wing and an analysis was 
conducted to determine the wing dimensions that provide optimum 
performance. 

The SPICA airfoil exhibits desirable characteristics for Reynolds 
numbers between 150,000 and 250,000. Thus, chord values ranging 
from 10 inches to 14 inches were selected. Given an estimated 
planform area of 9.5 square feet, based on minimum take-off 
requirements for a weight of 6.5 pounds and CLmax of 118, the range 
of corresponding spans was 6.97 feet to 13.68 feet. These 
dimensions resulted in a range of aspect ratios from 7 to 13 to be 
examined. The effects of varying aspect ratio on the drag polar, lift- 
to-drag ratio, and lift curve slope are illustrated and explained in 

Appendix A. The aerodynamic analysis revealed that a larger aspect 
ratio is desired to decrease induced drag, increase the maximum lift- 

to-drag ratio, and increase the lift coefficient for a given angle of 

attack. Minimizing drag decreases the amount of fuel required for 

flight, and therefore minimizes the operating expense. The 
maximum lift-to-drag ratio signifies that the aircraft is flying most 
efficiently and expending the least amount of fuel. A large lift 

coefficient at a small angle of attack is desired so that an adequate 
margin of safety exists between the angle of incidence and the stall 



angle. However, increasing the aspect ratio also increases the weight 
of the wing. 

Another important measure of merit is a minimum ratio of 
wing weight to lift. Minimizing the wing weight ensures that fuel is 
spent lifting the greatest percentage of cargo weight, thus 
maximizing the revenue earned each flight. Since the required lift 
of the wing is defined as the overall weight of the aircraft, our goal is 
to reduce the bulk of the wing while minimizing the induced drag 
that the wing develops at cruise conditions. However, these two 
objectives are contrary. 

As can be seen in the drag equation, Cd=Cdo + Cl 2 McAR (Ref 1), 
the wing’s effect on drag can only be diminished by increasing the 
aspect ratio and/or increasing the wing efficiency, since the parasite 
drag and lift coefficient are effectively defined by the aircraft’s total 
weight for a given a planform area and cruise speed. In order to 
achieve a maximum efficiency, the wing must be tapered, which 
would necessitate the construction of many different rib sizes and 
would significantly increase construction costs as stated previously. 
A more effective means of reducing drag is to increase the aspect 
ratio. However, since the effective lift force on higher aspect ratio 
wings is further from the fuselage, the bending moment at the root, 
and therefore the stress on the wing structure, is greater for higher 
aspect ratio wings. In order to accommodate this increased stress, 
the structure of the wing must be reinforced, making it heavier. As a 
simple model of this trade-off, the wing was modeled as the 
minimum rectangular spruce spar necessary to maintain a load factor 
of 2.5. Figure A.4 in Appendix A relates the approximation of the 
relation between wing weight and drag as the aspect ratio varies. 
This study determined that an aspect ratio of approximately 9.5 
provides the desired aerodynamic characteristics at an acceptable 
spar weight. 

A polyhedral wing configuration was designed to minimize the 
weight of the wing by utilizing continuous members across the 
fuselage and lighter joints near the tips. The polyhedral 
configuration functions like a dihedral wing by allowing the Prime 
Mover control authority without the additional weight and cost of 



ailerons. For ease of construction, the wing is divided into three 
sections with a 3.89 foot middle section and two 3 foot additions 
mounted at a 17 degree angle. This results in a 9.62 effective span 
and, with a one foot chord, a 9.62 effective planform area. Therefore, 
the effective aspect ratio is also 9.62. 

The two-dimensional lift curve slope for the airfoil was 
corrected for three-dimensional effects using the relation dCL/da = 
a D /( 1 +57.3a 0 /Tce AR) where a o =.0906/degree was determined from 
the lift curve slope for the two-dimensional airfoil (Ref 1). The 
efficiency factor was estimated to be .746 for an aspect ratio of 9.62. 
This value was determined by averaging the efficiency factors 
obtained from the two methods outlined in Appendix B. The lift 

curve slope for the wing is illustrated in Figure 4.3. Based on the 
given stall angle of the airfoil, the maximum lift coefficient was 
conservatively estimated to be 1.18. The lift coefficient required for 
cruise conditions was calculated to be .65 using the relation 

C L =2W/pVoo 2 S where the lift equals the weight of the airplane during 
steady, level flight (Ref 1). Therefore, the wing will be mounted at a 
7 degree angle of incidence to provide sufficient lift. An additional 
concern was tip stall. The change in angle of attack for the outboard 
wing sections may be approximated as Aa=BT (Ref 6). At maximum 
yaw angle of 10° and effective dihedral angle of 13.6°, one wing tip 
will at +2.3° while the other will be at -2.3° in reference to the center 

wing section. In a steady, level turn the maximum total angle of 

attack is 9.3°, an angle well below the stall angle. 

The variation of the lift-to-drag ratio is illustrated in Figure 
4.4. The maximum lift-to-drag ratio is 17.80 at an angle of attack of 
8.6 degrees. During cruise, the airplane flies at a lift-to-drag ratio 
only 2% less than the maximum value. Thus, the Prime Mover will 
fly efficiently at cruise conditions. 















4.3 Drag Prediction 

The drag for the entire aircraft was estimated using the 
relation for the drag coefficient Cp = Cdo + CL 2 MeAR where the first 
term is the parasite drag coefficient at zero lift and the second term 
includes both induced drag and the contribution of parasite drag due 
to lift (Ref 1). The method used for the estimation of Cdo was the 
drag breakdown method based on wetted surface area outlined in 
Jensen’s thesis. The total drag was found by adding the contributions 
from each component according to the formula Cdo= SCfS we t/S re f 
where a value of Cf=.0055 was obtained from empirical data and the 
wing planform area was used as S rc f (Ref 3). See Appendix A for 
detailed discussion of the drag breakdown. 

Given the values for aspect ratio, efficiency factor, and parasite 
drag coefficient the equation for the drag polar is Cd = .0179 + 
.0441 Cl 2 . Figure 4.5 shows the quadratic relation between the drag 
coefficient and the lift coefficient. 



Table 4.2 

Drag Component Breakdown 


Component 

CfSwct/Sref 

Wing 

.008685 

Fuselage 

.004156 

Horizontal Stabilizer 

.001616 

Vertical Stabilizer 

.000769 

Landing Gear 

.000310 

Interference 

Add 5% to Cdo 

Roughness and Proturberances 

Add 10% to C Do 

Total Cdo 

.0179 
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5.0 Propulsion System 

5.1 System Selection and Performance Predictions 

The Exodus Prime Mover requires a propulsion system that 
conforms to the following requirements: 

1. Range > 8600 feet. 

2. Velocity £ 30 feet per second. 

3. Take-off distance £ 60 feet. 

4. Cost as low as possible. 

To achieve these goals, the Prime Mover requires a system that works 
well together through all portions of the flight. The first step in 
assembling this system was to find a motor capable of providing enough 
power to enable the Prime Mover to meet the takeoff distance 
requirement. At the same time, the motor needed to be as small as 
possible to keep the weight of the system down. 

It was determined that fuel costs for the Prime Mover was 73% 
the total life cost of the aircraft. A study of current draw versus fleet 
life cost per volume moved showed that for every one amp increase in 
the current draw the fleet life cost per volume moved increased 15%. 
For these reasons, minimizing the fuel consumption became a prime 
concern for Group Exodus. The Prime Mover's, propulsion system 
consumes 86% of the total fuel per flight during cruise based upon 
estimates of required power (Figure 5.1). Therefore, it was decided that 
the studies involved in decreasing the fuel consumption would require 
an in-depth study of the cruise portion of the flight. 

Three motors were studied, the Astro Cobalt 05, 15, and 25. By 
looking at designs from previous years, it was noted that most aircraft 
were lighter than our aircraft and used either the 05 or the 15. It was 
therefore assumed that the 05 would be too small for our aircraft and 
was ruled out immediately. The 15 was analyzed to see if it would 
provide enough power for take-off. 



Figure 5.1 



The take-off program and the TK Solver program, electric motor 
performance, were used for the initial study. Because the motor needed 
to be ordered early in the design process, the study was done with 
many initial guesses for the required design parameters. This study 
showed that the 15 would provide enough power for take-off, even if 
the design parameters were to change. Cruise conditions were not a 
concern since there was more than enough power available from the 15. 
Therefore, the Astro Cobalt 15 was chosen because it provided the 
required power, weighed 30% less, and cost approximately 20% less 
than the Astro Cobalt 25. 

During the take-off portion of the flight, the motor will be running 
at 15,192 revolutions per minute at an efficiency of .73. During climb 
the motor is running at 16,195 revolutions per minute at an efficiency 
of .75. Finally, the motor will work for the rest of the flight at the 
cruise condition of 11,186 revolutions per minute at an efficiency of .55. 
These figures were based on manufacturers specifications and computer 
analayses found in Appendix G. A graph of the motor efficiencies 
versus the rpms can be seen in figure 5.2 below. 
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Figure 5.2 

Motor Efficiency versus RPM 
Astro 15 



The next piece of hardware that needed to be ordered was the 
batteries. It was determined from the initial studies on the Astro 15 

that during take-off the engine would require 14.4 volts to provide the 
power required. This voltage requires 12 batteries connected in series, 
since each battery had an individual voltage of 1.2 volts. 

The next step was to determine the capacity of the batteries 
needed for the Prime Mover. The objective of Group Exodus was to 
have a range larger than 8600 feet, which includes diversion and loiter. 
For take-off and climb the motor uses a total of only 20 milliamp hours. 
At 30 feet per second the current draw for the motor was 4.45 amps. 
To cover the range of 8600 feet the motor would use a total of 354 
milliamp hours, for a total of 376 milliamp hours used for the entire 
flight. 

Since each aircraft would be making two flights a day, it was 
decided to use batteries with the capacity to make this round trip flight 
without needing to change the batteries. This would provide for a 
decrease in operational costs as well as a decrease in the length of time 
on the ground. Therefore, a battery capacity of 752 milliamp hours was 
needed. The nearest capacity battery to this required amount was a 

1000 milliamp hour battery, which provides more than enough fuel for 
the Prime Mover. 




5.2 Propeller Design 

Determination of a propeller for the system proved to be a 
difficult task. Picking a propeller to keep the current draw during 
cruise as low as possible was the main task, for the reasons stated 
previously. The fuel consumption and total efficiency of the system was 
studied during this portion of the flight for different propellers. 
Unfortunately, it was determined that an efficient propeller for the 

cruise portion of the flight didn't necessarily allow the aircraft to take- 
off in the required distance. Therefore, a study was made of several 
different propellers to determine the best for the entire flight regime. 

A major tool used in this selection was the propeller program 

written by Barry N. Young. This program predicts various performance 

characteristics, such as coefficient of thrust, power, and the efficiency, 

for a particular propeller design. This is done by use of the blade 
element theory which included considerations for induced velocity and 
tip losses, as well as lift and drag coefficient adjustments for the 
working Mach and Reynold's numbers. The use of this program 
requires the user to input specific data on the propeller, namely airfoil 
sections, chord, thickness, and angle of the blade at different radial 
positions. An example of the propeller input data and the output from 
the program can be seen in Appendix D. It is important to note that the 
airfoil section selected for this program should be the NACA44XXLOWRE 
because it is a low Reynold's number airfoil. The reason for this is that 
the propeller is working at very low Reynold's numbers, approximately 
60,000 at 75% of the diameter, which can be seen on the ZingerJ 11-5 
performance estimate sheet in Appendix D. It is believed that by doing 
this the performance predictions for the propeller will be more 
accurate. 

This program was run for several propellers ranging in diameter 
from 8 inches to 12 inches and in pitch from 4 inches to 6 inches. The 
data was compiled onto Cricket Graph in order to obtain more 
information to be used in the TK Solver program, electric motor 
performance. An example of the input data and the rules used for the 
TK Solver program can be seen in Appendix E. This program requires a 
value for the coefficient of thrust and power for an advance ratio of 
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zero. This value is not available from the propeller program therefore 
the values had to be extrapolated for an advance ratio of zero. Once 
these values were found the electric motor performance was utilized to 
find the efficiencies and the current draw during cruise for the studied 
propellers. The results can be seen in Figures 5.3 and 5.4. 

From this it can be seen that the best choices would have been 
either the 8-6 or the 10-4. Unfortunately, neither of these propellers 
would enable the aircraft to takeoff in the required 60 feet. Therefore, 
the ZingerJ 11-5 was chosen for our aircraft. 

During take-off the propeller will work at an advance ratio of .2 
and an efficiency of .61. During climb the propeller will increase it's 
advance ratio to the cruise condition of .42 and work at it s maximum 
efficiency of .745. This can be seen in the propeller efficiency versus 
advance ratio curve, figure 5.5. 


current 
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Figure 5.3 

Current Required at Cruise 

Weight=6.5 lbs. 
Velocity=30 ft/s 
Aspect Ratio=9.5 
Wing Area=9.5 ft A 2 
Prop Diameters 11 in 


8-6 9-6 10-4 10-6 11-5 12-4 12-6 

propeller 



w 


&» 


Figure 5.4 

Efficiency at Cruise 



8-6 9-6 10-4 10-6 11-5 12-4 12-6 

propeller 


Efficiency 


Figure 5.5 

Propeller Efficiency versus Advance Ratio 

ZingerJ 11-5 
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5.3 Engine Control 

The motor will be controlled by the Tekin speed controller. 
Throughout the flight the Prime Mover will be operating at two 
different throttle settings. At take-off the Prime Mover will be at full 
throttle(14.4 volts). This results in a maximum rate of climb of 5.4 feet 
per second for the aircraft at a forward velocity of 30 feet per second, 
the Prime Mover's designed cruise velocity. Once the cruise altitude is 
obtained the Prime Mover will be throttled back to 67% full throttle(9.6 
volts). This information is shown on the power required and power 
available versus flight speed curve in Figure 5.6. 





6.0 Preliminary Weight Estimation Detail 


6.1 Component Weight Estimates 

Based upon an economic analysis of the available Aeroworld 
overnight delivery market. Group Exodus set out to design a cargo 
aircraft capable of transporting a 2 pound payload at a velocity of 27 
ft/s. A rough estimate of the airplane's total weight, 6 pounds, was 

obtained by extrapolating from last year's AE441 large-scale, battery- 
operated aircraft concepts, which were designed for smaller payloads 
and weighed an average of 4.6 pounds. Using this initial weight 

estimate, we approximated the weights of the various components to 
verify the initial guess. 

Based upon the total weight, it was determined that a wing 
planform area of 9 ft 2 would be needed to maintain steady, level flight 
at cruise conditions with an estimated Cl of .8. From this planform area, 
an estimate of the wing weight was made based upon a linear 

extrapolation of the densities of previous wings. As is shown in Table 

6.1, the densities of previous wings (wing weight/planform area) vary 

between .01 and .026 oz/in 2 . Since the chord was estimated at 1 foot, 

we felt that the wing of the Behemoth Apteryx, which was slightly 

overdesigned, would be most like our wing design. With a few 

technological advances, we felt that we could design a wing of roughly 
twice the span at the same density, .011 oz/in 2 , for a total wing weight 
of about 14.5 ounces. 
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Table 6.1 

Wing Density of Existing Aircraft 

Aircraft °S (so in) | Weioht (oz) [Density (oz/in2J_ 

.011 
.017 
.012 
.010 
.026 
.014 

The weight estimate of the fuselage was determined in a similar 
manner as the wing weight. Since it was desired that the interior of the 
cargo bay to house a line of 4 inch boxes, the exterior dimensions of the 
fuselage were estimated at 4.5" x 4.5" x 60". The fuselage of last year's 
Initial Guess had about the same cross-section, and weighed 11.8 oz. for 
a 51 inch length. Scaling its weight for our 60 inch fuselage, an 
approximate fuselage weight of 14 ounces was obtained. 

In choosing a propulsion system, it was noticed from the engine 
specifications that there is a large weight penalty when upgrading from 
an Astro 15 to an Astro 25; therefore, it was an objective to stay with 
the lighter Cobalt 15 engine. For the initial weight estimate of the 

Prime Mover, the avionics and propulsion system weight estimates 

provided in the spec sheet were used. 

Finally, since the empennage and landing gear contribute so little to 

the overall aircraft weight, a very rough estimate on the weights of 
these components was made based upon comparable components 
utilized in previous aircraft. These weights, along with the weight 

estimates and c.g. locations of all of the other components of the Prime 
Mover, are listed is Table 6.2a. From these approximate component 

weights, the total weight estimate of the Prime Mover came to about 6.8 
lbs. Since this figure is slightly higher than the initial estimate, the 
wing area was increased to about 9.5 ft* with a weight of 15 oz. These 
initial weight and location estimates are very close to the actual figures 
for the technology demonstrator (Table 6.2b), which weighed only 6.7 
pounds when fully loaded. 


840 

1010 

1000 

1000 

1440 

842 



Table 6.2a 

Weight Estimates and C.G. Locations of Aircraft 
Components 


Component 


Weight (oz) 


x position (in) 


y position (in) 


z position (in) 


fuselage 

engine 

empennage 

Landing Gear (front) 

Landing Gear (rear) 

Cargo 

receiver 

servos 

system battery 
speed controller 
wing (left) 
wing (right) 
battery 1 

battery 2 

I total weight 


1 4 

-30 

10.25 

- 2 

5 

-57 

5 

- 3.5 

2 

-58 

32 

- 30.75 

0.95 

- 4.95 

1 .2 

- 4.95 

2 

- 4.95 

1 .8 

- 4.95 

7.5 

- 32.25 

7.5 

- 32.25 

10.02 

- 49.5 

10.02 

- 49.5 

109.24 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

30 

-30 

1 

- 1 


0.1 

4 . 


- 0 . 

1 .25 
-1 .5 
0.5 
6 
6 

2.5 

2.5 


o to h*- co o in 



Table 6.2b 

Final Weights and C.G. Locations of Aircraft Components 


Component 

Weiaht (oz) 

x Dosition (in) 

v Dosition (in) 

z position (in) 

fuselage* 

16.9 

-31 .8 

0 

0.1 

engine** 

10.25 

- 3 

0 

0 

empennage 

2.8 

-58.5 

0 

4.5 

Landing Gear (front) 

4.6 

-19 

0 

- 5 

Landing Gear (rear) 

1 

1 -59 

0 

-4.5 

Cargo 

32 

-32 

0 

0 

control package*** 

6.6 

- 9 

0 

0 

wing (left) 

4.6 

-32 

38 

7.6 

wing (mid) 

6.2 

-32 

0 

3.3 

wing (right) 

4.6 

-32 

- 3 8 

7.6 

battery 1 

9 

-51 

1 

3 

battery 2 

9 

-51 

- 1 

3 

1 total weight 

107.55 





* includes two control rods and 4 ft battery cable (2 oz total) 

** includes geared engine and propeller 

*** includes receiver, servos, system battery, and speed controller 



6.2 Center of Gravity Location and Travel 


The center of gravity of the Prime Mover is located 1.2 inches above 
the centerline of the fuselage and 30.5 inches behind the propeller 
(Figure 6.1). Due to the placement of the batteries atop the rear 
fuselage, the center of gravity of the Prime Mover is practically 
independent of cargo weight. 


Figure 6.1 

Position of the Center of Gravity 



\ 


1 .2" from 
centerline 


As shown in Figure 6.2, for the range of possible cargo densities, the 
center of gravity of the entire aircraft deviates only 1% of one chord 
length. Since the Prime Mover possesses a continuous cargo bay, the 
C.G. of the cargo can be held effectively constant by judicious 
placement of cargo. The center of gravity may move from design 
conditions more significantly, however, due to errors in component 
weight estimates; yet, as can be seen in Figure 6.2, this off-design 
performance accounts for a C.G. movement of at most 8% of one chord 
length. 


C.G. position from wing L£. (in index units) 
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Figure 6.2 
Off-Design Performance 



Airplane weight (lbs) 


Design conditions 
-30% Design Weight Error 
+30% Design Weight Error 



7.0 Stability and Control 


Stability can be divided into three major areas of concern: 
longitudinal, directional, and roll stability. 

7.1 Longitudinal Stability 

For an aircraft to have longitudinal static stability, the pitching 
moment curve must have a negative slope (ie Cma < 0). Appendix F.l- 
F.3 shows the contribution of each component and the assumptions 
involved. The neutral point of the aircraft is where Cma = 0. It is 

defined by: 

Xnp/c = Xac/c - Cmaf/Claw + T|h*Vh*Clah/Claw*(l-d£/da) 

The static margin is a measure of the degree of stability and is defined 
by: 

Static Margin = Xnp/c - Xcg/c 

The larger the static margin, the more stable an aircraft will be. 
Through an analysis of similar aircraft from previous years, it was 
determined that a static margin of approximately 15% should allow safe 
flight without requiring large control inputs for maneuvering. The 
static margin values for typical aircraft are between 5° and 10°. The 
reason for the higher values for this type of aircraft is the larger 
response time needed by a remote pilot. 

There are many factors that effect longitudinal stability, but only 
a few which are not set by other requirements. These are: 

1. Center of Gravity (CG) location 

2. Wing Position 

3. Horizontal stabilizer incidence angle 

4. Horizontal stabilizer moment arm 

5. Horizontal stabilizer area/aspect ratio 
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The CG location plays a major role in longitudinal stability. The 
more the CG travels, the more the handling characteristics of the Prime 
Mover will vary. In order to reduce this variation, the battery CG was 
used to make the Prime Mover CG approximately independent of 
payload. Therefore, the CG location, though not set by requirements 
outside of stability and control, was set to limit the variation in handling 
characteristics of the Prime Mover. The amount of CG travel between 
full and empty payload cases is only one tenth of an inch. A graph of 
Static margin versus CG for the Prime Mover's forward and aft CG limits 
is shown in figure 7.1 (please note that inches from the leading edge 
should read inches aft of the leading edge). It can be seen that static 
margin is strongly influenced by CG location. The arrows on this and all 
figures in this section represent the present design values. 

FIGURE 7.1 

Effect of CG Travel on Static Margin 



Xcg (inches from leading edge) 
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As mentioned earlier, the angle of attack of the fuselage reference 
line and horizontal stabilizer incidence angle are both desired to be zero 
at cruise in order to reduce drag. Assuming this, for the aircraft to be in 

equilibrium during cruise, the wing lift force and pitching moment 

should be balanced to result in zero moment about the CG. Altering the 
wing position is very useful in controlling this condition. Appendix F.4 
shows how the wing position was set. The CG position was determine 
as described above and the coefficient of lift and chord were set by 
cruise conditions. However, there was no data for Cm ac for the SPICA 
airfoil. It was observed that symmetric airfoils have Cm ac values of 
zero and cambered airfoils have Cm ac values proportionate to the 
amount of camber. It was also observed that the more camber, the 
more negative the value of the oc c i = o on the lift curve slope. Combining 
this information, it was assumed that the SPICA would have a similar 
Cm ac as an airfoil that has the same Ct c i=o. Therefore, Cm ac was 
assumed to be -.05. The resulting wing position placed the aerodynamic 
center at 30% chord point. This is in agreement with traditional 
conventional aircraft. 

A higher angle of incidence of the horizontal stabilizer creates 
more drag at cruise cruise conditions. To avoid this unnecessary drag, it 
was desired to have the horizontal stabilizer attached at zero angle of 
incidence. This would also make construction simpler, hence reduce 
the construction hours and cost. However it was not possible to have 

the tail at zero angle of incidence. The angle of incidence for the 

horizontal stabilizer was set at +5°. The reason for the large angle is 
directly related to the small moment arm which will be discussed next. 

As the horizontal stabilizer moment arm increases, the 
effectiveness of the horizontal stabilizer increases, hence the smaller the 
horizontal stabilizer area required to have the same effect. The smaller 
horizontal stabilizer reduces both drag and weight. Therefore, the 

horizontal stabilizer was placed as far aft as possible on the fuselage. 
Since the CG is so far aft on the fuselage, the moment arm is small 
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relative to the size of the fuselage. This will be a major hindrance in the 
design process. 

The stability of an aircraft can be increased by increasing the 
horizontal stabilizer area; however, if the horizontal stabilizer is too 
large it will require a large elevator or elevator deflection to achieve 
CLmax. Flying at Clmax allows the aircraft to fly at the slowest possible 
speed (V s tall) making landing easier. Horizontal stabilizers with higher 
aspect ratios are more effective. Large bending moments associated 
with higher aspect ratios arise because the horizontal tail is a simple flat 
plate and offers only small amounts of stiffness about the roll axis. The 
sensitivity of Static Margin to horizontal stabilizer area and aspect ratio 
is shown in figure 7.2. The Prime Mover has a horizontal tail area of 
1..44 ft^ and an aspect ratio of 3.25. This results in a static margin of 
15.6%. 


FIGURE 7.2 

Effect of Horizontal Tail Area and 
Aspect Ratio on Static Margin 
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The elevator was sized to allow achievement of CLmax. The 
deflection angle was limited to +/- 15°. This is because large deflection 
angles result in high drag. Using an elevator size equal to 18.75% of the 
horizontal tail area, a graph of Cm versus a for maximum elevator 
deflection is shown in figure 7.3. Since the CG travel is basically zero, 
there are no real forward and aft limits for the CG. As long as the Cargo 
Cg is placed at the center of the cargo bay, the overall CG of the Prime 
Mover will be the same. From this it can be seen that an elevator 
deflection of 12° will provide the Prime Mover with Clmax. Figure 7.1 
can be used to determine the allowable limits on CG motion. The limits 
are set by the neutrally stable condition (SM = 0%) and a static margin 
25%. It shows a 2.7 inch CG travel, however, the Prime Mover will want 
to maintain a tolerance of 0.5 inch to maintain adequate handling 
qualities. This represents the stupidity tolerance of the Prime Mover, 
since it will only be of concern if the cargo is loaded improperly. The 
elevator will be controlled by a servo in the avionics compartment. The 
control cables will run along the sides of the cargo bay. 




7.2 Directional Stability 

For an aircraft to have directional static stability, it must have a 
positive yawing moment curve slope (ie CnB > 0). The main contributors 
to CnB are the fuselage and the vertical stabilizer. Appendix F.5-F.6 
show how each component's contribution was determined. 

The important parameters for directional stability that were not 
set by other considerations are: 

1. vertical stabilizer moment arm 

2. vertical stabilizer area/aspect ratio 

Vertical stabilizer moment arm was desired to be as large as 
possible to reduce the size of the required vertical stabilizer, thus 
reducing drag and weight. Therefore, it was placed as far aft as 
possible. Again, the relatively short moment arm of the Prime Mover 
adds difficulty to the design. 

Large vertical stabilizer area increases weathercock stability. 
However, if the vertical stabilizer is too large it will have a problem 

with large control input requirements analogous to the large horizontal 

stabilizer. The vertical stabilizer aspect ratio increases the effect of the 
vertical stabilizer, but is limited by structural concerns. 

7.3 Roll Stability 

For an aircraft to have roll static stability, it must have a negative 
banking moment curve slope (ie C1B < 0). The main contributor to C1B is 
effective dihedral angle (EDA). C1B was calculated by the method shown 
in Appendix F.7. 

Since the Prime Mover is employing a rudder/dihedral 

combination, as opposed to ailerons, to meet the turning requirement, it 

will require a large EDA. Therefore, EDA selection was driven by the 
turning requirement instead of stability. At the design cruise velocity, 
and assuming a 25°/second roll rate, the Prime Mover requires 13.2° 
EDA, which results in a 17° panel dihedral since the Prime Mover's 
break point is at 40% of the span (appendix F.8). Using this information, 
the Prime Mover would require a 10° yaw angle to meet the turning 
requirement (appendix F.9). The rudder was sized to be able to achieve 
a 10° yaw angle. The deflection angle was again limited to +/- 15° for 
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the same reasons. A graph of Cn versus B (figure 7.5) for a rudder 

area of 56.7% of the vertical tail area shows that a rudder deflection of 

15° will satisfy the turning requirement. The rudder will also be 

controlled by a servo located the the avionics compartment, with a cable 

running along the sides of the cargo bay. 


FIGURE 7.5 

Yawing Moment versus Yaw Angle 
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8.0 Performance Estimations 

8.1 Take-Off and Landing Estimates 

Estimation of take-off performance was a complicated yet crucial 
task that became the determining factor in the sizing of the aircraft 

components. The take-off analysis began with the use of a computer 
code, listed in Appendix G. Following the selection of the motor and 
batteries, it was determined that the most influential specification on 
take-off distance was the propeller type. To insure realistic values of 
the propeller coefficients, corrections were made for velocity and tip 
losses, Reynold’s number, and airfoil section. Because of the mission-set 
objective, requiring the Prime Mover to complete its take-off run under 
60 feet, the propeller was limited to a minimum diameter of 1 1 inches. 
The final propeller chosen was the ZingerJ 11-5. 

Once the sizing of the aircraft was complete, cargo weight and 

coefficient of friction were varied to test the sensitivity of take-off 

distance. Figure 8.1 shows the variation of take-off distance with 

weight and coefficient of friction. The Prime Mover was designed to 
take-off under the required distance, carrying the maximum possible 
cargo weight, and subject to the estimated coefficient of friction of the 
targeted run-ways (.15). For larger values of coefficient of friction the 
cargo weight need be reduced. For a more detailed explanation of the 
take-off calculations, see Appendix G. 
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Figure 8.1: Variation in Take-Off Distance 

with Cargo Weight and Coefficient Friction 



Landing distance was another important estimation. Without the 
use of brakes these distances can be quite large. Fortunately, the value 
for coefficient of friction for the run-ways was also large, helping to 
reduce the landing distance. Estimation of ground roll distance during 
landing was made with the use of the equations found in Reference 5. 
Exact calculations are shown in Appendix G. The landing distance for 
the Prime Mover was found to be 44 feet at maximum cargo weight, 
enabling landing in all targeted cities. Although no braking system will 
be employed during the flight test, plans for construction of brakes 
should be completed for manufacturing of the fleet to allow for changes 
in coefficient of friction. 

8.2 Range and Endurance 

Range and endurance estimates were found using the computer 
code listed in Appendix G. Figure 8.2 shows the values for range and 
endurance of the fully loaded Prime Mover varying with velocity. The 
final results far exceeded the mission objective of a minimum range of 
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8900 feet. The reason for this occurrence was due to the 1000 
milliamp-hour capacity of the selected batteries. A much lower battery 
capacitance should have been used in order to reduce the amount of 
over-shoot. 


Figure 8.2: Range and Endurance 
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Under the selected cruise velocity of 30 ft/s, the range and endurance 
were calculated as about 24,000 feet and 13.2 minutes, respectively. 
The maximum range occurred around 31,000 feet at 45 ft/s with a 
corresponding endurance of 11.5 minutes. At 33 ft/s, the maximum 
endurance of 13.5 minutes was discovered with a range of 27,000 feet. 
Another important consideration was the variation in range with 
payload weight. Figure 8.3 shows this trend with values of cargo weight 
from empty to maximum. Again, with the 1000 milliamp-hour 
batteries, there was no concern with the inability to meet the design 
objectives. 
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8.3 Power Available and Required Summaries 

The power available calculations were completed for various 
voltage settings, using the computer code in Appendix G. Results of this 
study are shown in Figure 8.4 for the fully loaded Prime Mover. The 
significant features to note of these results are the values of minimum 
and maximum velocity and the voltage setting for the desired cruise 
velocity. Maximum and minimum velocities occur where the power 
available for maximum voltage setting intersects the power required 
curve. These values were found to be about 10 ft/s and 51 ft/s, 
respectively. 


Figure 8.4: Power Available and Required 

for varying voltage settings 



At the cruise velocity of 30 ft/s, the voltage setting was about 9.5 volts. 
Also shown in the figure is a power available curve for one other 
voltage setting. This example shows how to determine the voltage 
setting for various desired cruise velocities and rates of climb. 
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8.4 Climbing and Gliding Performance 

As a result of the power available and required study the rate of 
climb for various velocities and voltage settings was determined. Of 
particular importance was the maximum rate of climb, occurring at the 
maximum voltage setting, found to be 5.4 ft/s. This rate of climb occurs 
at the desired cruise velocity of 30 ft/s. The ratio of horizontal to 
vertical distance traveled during maximum rate of climb is 5.56, 
meaning a significant amount of ground distance must be covered to 
increase in altitude. To illustrate this fact, assuming the maximum rate 
of climb can be achieved directly after take-off, about 140 feet will be 
needed to climb to a 25 foot height, over twice the distance needed for 
take-off ground roll. 

In case of such emergencies as engine failure, it was important to 
investigate the Prime Mover's gliding performance. The minimum glide 
angle, giving the maximum horizontal distance covered in gliding flight, 
occurs when the lift-to-drag ratio is greatest. The minimum glide angle 
was found to be 3.33°, giving a horizontal distance travelled of 430 feet 
from a starting altitude of 25 feet. 

8.5 Catapult Performance Estimate 

The series of flight tests includes a catapult test to test such 
characteristics as aerodynamics and stability. To estimate the predicted 
performance during this test, a catapult program, written by Kevin 
Costello, has been utilized. 

One particular problem encountered during the calculations was 
the difficultly in obtaining a relatively large range while touching down 
approximately parallel to the ground. Many of the initial flight 

conditions tested either resulted in the inability to take-off or the 
occurrence of a nose dive. Employing an elevator deflection helped to 
alleviate this problem. Figure 8.5 presents the results of one particular 
catapult estimation in which the variable parameters were as follows: 
Distance between hard points : 20 feet 
Catapult deformation : 30 feet 
Elevator deflection : -0.15 rad 
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Figure 8.5: Catapult Estimation 



Horizontal Distance Covered 

The catapult range measured from the pins was 35 feet with the angle 
of the fuselage reference line at touch down at 0.009 rad. The 
maximum altitude achieved was about 2 feet. 
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9.0 Structural Design Detail 

9.1 V-n Diagram 

The Prime Mover is not only designed to withstand the expected 
flight conditions, but also, to handle the unexpected loading 
environments of Aero World. Figure 9.1 shows the load limitations of 
the aircraft versus the aircraft velocity for an estimated weight of 
6.5 lb. The maximum load capacity during flight maneuvers is 2.5 g, 
and during landing is 4.0 g. The maximum velocity of the Primer 

Mover is 51.0 ft/sec. However, due to the Aero World constraints on 
noise, the Prime Mover will not fly faster than the speed of sound, 30 

ft/sec. 

9.2 Basic structural components, substructures, and assembly 

The elimination of the need for ailerons makes a dihedral wing 

configuration very attractive due to savings in both construction and 
maintenance costs. However, the bulky reinforcement at the root of the 
dihedral wing means a possible weight penalty. Use of a polyhedral 
wing eliminates the middle joint, allowing a continuous main spar to 

carry bending stress. Since the wing joints are further out on the wing 

and since the bending moment decreases as the square of the distance 
from the root, the polyhedral wing utilized on the Prime Mover 

provides a lighter solution for aileron elimination than does a dihedral 
configuration. 
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The major stresses in the wing are carried through the C-beam 
main spar. The spruce spar caps, in particular, are responsible for the 
integrity of the wing under the stress caused by lift forces. For these 
caps, spruce, due to a significantly higher yield stress than balsa, was 
found to weigh less than the amount of balsa necessary to prevent 
failure at the root of the wing. Balsa spar webs will also be included in 
the middle-section of the wing, where bending moments are the 
highest, to prevent excessive shear stresses on the ribs. The leading 
edge and trailing edge spars are of balsa construction and carry very 
little structural load. The ribs are also made of balsa and primarily 
maintain the integrity of the airfoil shape. Since the monokote is most 
likely to droop at the leading edge, where there is a high pressure 
gradient, half ribs are employed every six inches to increase the 
aerodynamic effectiveness of the wing. Toward the trailing edge, less 
structural mass is wasted supporting level sections of skin, which 
accounts for a rib weight savings of about 10% over a more closely 
spaced arrangement of full ribs. Wherever possible, the ribs are cored 
out to further reduce their mass. 

Figure 9.2 


Middle section of wing structure 



The most complicated piece of the wing is the joint between the 
main wing section and each canted, exterior wing section. These 
exterior sections must be detachable so that the Prime Mover will meet 



storage requirements; yet, the joints must also be strong enough to 
weather a hard landing. Since an angled joint must be employed to 

connect the main spar of each section, a light plywood must be used, at 
a small weight penalty, to avoid failure due to shear along the grain. 
These plywood joints are affixed to the exterior sections and fit into 
sleeves between the spar caps of the middle wing section. Pins at the 
leading and trailing edges and an additional rib along the joint ensure 
that the wing sections are flush with each other, eliminating lift 
degradation due discontinuities in the skin. 

At the center of the middle wing section, a reinforced spar web 
drops down and anchors the wing to the top of the fuselage. Screws at 
the trailing edge of the ribs which flank the fuselage anchor the wing to 
the aircraft at the proper incidence angle. 

The size of the empennage is determined to optimize the control 
and stability of the aircraft. (See chapter 7). The empennage is 
designed to be easily removable and light weight in order to meet 
storage requirements and weight objectives. 


Figure 9.3 
Empennage Design 



With the exception of the mounting hooks, the entire empennage 
is made of balsa. The size of the balsa wood was chosen to ensure 
handling strength while minimizing the weight. The configuration of 
the empennage components are shown in Figure 9.3. The empennage 
will be connected by screws in the aft section and a metal mounting 
hook in front. 
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The fuselage was designed to keep a minimum frontal area while 
allowing room for either the two or four inch cubes in the cargo space as 
illustrated by Figure 9.4. The circular holes in the upper corners 
represent the location of the control cables with respect to the cargo 

space. 


Figure 9.4 

Fuselage Cargo Spacing 


4.625" 


4.275" 



The side of the Prime Mover fuselage structure is supported by a 
sawtooth truss design (Figure 9.5). This was chosen over right angle 
and left angle truss designs because it is stronger than either in both 
tension and compression as evaluated by the Exodus internal stress 
analysis program (See Appendix H). The location of the vertical cross 
beams were chosen to minimize the fuselage weight while withstanding 
4.0 g landing loads, 2.5 g cruise loads, and catapult loads of 2 g's. 
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In order to simplify construction, the beams where then moved to 
even intervals and retested. The amount of weight added by this step 

was only a fraction of an ounce. 

Figure 9.5 

Fuselage Truss Design 

The forward landing gear is located directly under the instrument 
package where the fuselage is already strengthened. Unfortunately, 

this is located far forward of the center of gravity, possibly producing 
an undesirable attitude if landing when the aircraft is subject to large 

yaw moments. To counter this situation the landing gear is given a 

wide base. 

9.3 Internal Configuration 

There are four major internal components. These are: 

1 . motor 

2. batteries 

3. avionics 

4. cargo 

The location of each component plays an important role in the 
performance of the aircraft. The internal configuration can be seen in 
figures 9.6 and 9.7. 
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Figure 9.6 



30.5" from prop 

Figure 9.7 

Instrument Package Configuration 



Motor 

The motor is located in the front of the Exodus to reduce blockage. 
This location also reduces the risk of ground contact since the Exodus 
has tail dragger landing gear. The motor is centered and slightly high 
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on the fuselage in order to align the motor with the CG of the Exodus 
and minimize the effects of the propulsion system on the moment about 
the CG. Another benefit of having the motor slightly high on the 
fuselage is that it reduces the required landing gear height. 

Batteries 

Since the fuselage, motor, and cargo (full payload) CGs are fixed, 
the batteries are the only remaining component with sufficient weight 
to significantly effect the total CG of the aircraft. It is desirable for the 
CG of the Exodus to be independent of the payload it is carrying so that 
the handling qualities are consistent. In order to achieve this it was 
necessary to make the CG of the complete aircraft with no payload 
match the CG of the cargo only. This required the batteries to be as far 
aft on the fuselage as possible. 

The batteries are on top of the fuselage in a separate 
compartment for several reasons. It provides easy access, thus lower 
maintenance cost. It also maximize internal volume. Another 

important reason is to provide cooling. The batteries get very hot and 
ventilation will be much better in a separate compartment where vents 
can be added. The wires from the batteries will run inside the main 
fuselage along the sides of the cargo bay. 

Avionics 

The Avionics is composed of five components: 

1. Receiver 

2. Speed Controller 

3. System Battery 

4. Servo #1 

5. Servo #2 

These components are located close to each other for easy access and to 
reduce the wire length between them. Easy access is a requirement 
because the entire system must be able to be installed in less than 30 
minutes. Each component will be on a separate platform, and each 
platform will be connected. This system will have a handle to provide 
easy removal of the total system. The servos will be connected to 
control cables for the elevator and rudder. These cables will run along 
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the sides of the fuselage in a similar manner as the battery wire. This 
requires the servos to be on top to allow for cable connection after the 
components are in place. The avionics are located forward on the 
Exodus to limit the load on the compartment since the aircraft will be at 
an angle of attack while on the ground. 

Cargo 

The Exodus has a single uninterrupted cargo bay to provide 
complete control capability of the cargo CG. The boxes will be supported 
by strips of wood running the length of the cargo bay. These strips will 
have holes every inch to allow for a peg to be inserted, thus 
constraining the boxes during non-full cases. The door to the cargo 

bay will have 1/8 inch guides to insure the cargo does not interfere 
with the control cables or battery wires running along the side of the 
compartment. 



10.0 


Construction Plans 


10.1 Major Assemblies 

As with any type of construction, difficulties can be anticipated 
and overcome by early planning. For this reason. Group Exodus has 
chosen a number of major assemblies which need extra attention. 

The first is the polyhedral joints. These joints will be connected in 
a sleeve joint assembly. The angle will be constructed into the joint and 
held by additional slanted rib supports. The joint will be made out of a 
light plywood which will slide between the spar webs and spar caps of 
middle wing section (Figure 10.1). The plywood will be designed to 
maintain a tight fit in the .25" by .9" webb area of the inner and the 
outer wing sections. 


Figure 10.1 
Polyhedral Joint 



The second obstacle is the wing-fuselage attachment. The wing 
will be attached to the fuselage by two means. Two hooks originating 
from the wing spar will attach to a strengthened fuselage cross beam 
near the outside of the fuselage. The back of the wing will be attached 
by screwing the interior ribs to the fuselage. The locations of the rear 
attachment are selected so as to produce the proper wing incidence 
angle. 

The third challenge is the empennage mounting which will be 
done in a fashion similar to that of the wing. The vertical tail hooks 
around the horizontal stabilizer which in turn is hooked to the fuselage 
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near the outer edges of the fuselage. This configuration minimizes the 
forces exerted on the fuselage support beam. The entire empennage is 
then fastened to the fuselage at the rear of horizontal and vertical tails 
just prior to the elevator and rudder by use of a large plastic screw 
(Figure 9.3). 

The instrument package is placed between the engine firewall and 
the instrument protection wall. Sufficient room is allowed for easy 
removal and reconnection of the instruments as necessary. The 
forward landing gear is connected directly under the instrument 
package. 

Figure 10.2 

Removable Instrument Package 



Connecting the rear landing gear offers some difficulty. An 
attempt must be made to give the rear landing rotational freedom in 
order to assist in ground handling. This freedom can be assured by 
attaching the landing gear in a pivot block made of plywood. Washers 
fastened to the gear will stop landing gear vertical motion while 
allowing turning capability. See Figure 10.4 
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Figure 10.3 

Rear Landing Gear Attachment 

Support Washer 



10.2 Complete Parts Count 

The parts count is based on the initial blueprint of the Prime 
Mover. Alterations will need to be added as construction begins. For 
this reason a price inflation of 10% is added to each component to meet 
unexpected changes in the prototype construction. 

Figure 10.4 is a listing of the initial part count. 

10.3 Assembly Sequence 

Assembly of the vertical tail, horizontal tail, wing sections, and 
empennage side trusses will all be done independently. At this point 
the aircraft's center of gravity can be verified. The final stages of 
construction demand extra attention to assure proper connection with 
the wing and empennage as well as proper placement of the battery 
and wing. The complete fuselage will be put together at this point after 
verifying all connection designs. 

The over all production time should be enormously decreased by 
this modular construction plan. With smart planning, production of the 
Prime Mover should easily be achieved in less then the one hundred 
and twenty man hour goal. 
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11.0 Enviromental Impact and Safety Issues 

11.1 Disposal of Each Component 

The revolutionary success of the Exodus Prime Mover will ensure 
that the first few hundred production units will be sent to 
miscellaneous museums throughout the world. 

Keeping environmental considerations in mind. Exodus plans on 
reusing or recycling all products to the end of their expected life cycle. 
For instance, all propulsion units are rechargeable and aircraft control 
equipment may be reused in other airplanes. 

Difficulties arise in recycling wood because all connections 
between members are joined by toxic glue. All salvageable wood may 
be recycled for any use (i.e. gardening mulch or firewood). Infected 

wood and plastic monokote materials must be handled separately as 
potentially hazardous materials. 

11.2 Noise Characteristics 

Unfortunately, Aeroworld is free of real world effects associated 
with the. effects of exceeding Mach one. The aircraft itself will be 

limited to a velocity of mach one; however, the propellor can not meet 
that limitation. The rotational speed necessary to propel any aircraft 

results in a tip speed far greater than 30 ft/sec. The noise elements 
related to the propellor speed can not be avoided. 

The actual hub station should be placed outside the residential 

area of city J in order to reduce the number of complaints related to late 
night noise. Other suggestions for noise reduction include detoured 
flight patterns as well as engine muffling or resizing. 

11.3 Waste and Toxic Materials 

Disposal of waste and toxic materials is generally subcontracted to 
waste disposal companies. In the event that these companies can not 
deal with glue or monokote materials City O, the New Jersey of 

Aeroworld, would be an ideal dump site. 
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12.0 Economic Analysis 

Probably the most important goal of this project was to maximize 
the return on investment. Consequently, the economic analysis was one 
of the most crucial studies. In fact, the economic strategy affected the 
mission and the final design more than originally expected. 

Before performing the economic analysis. Group Exodus adopted 
an economic strategy. We realized that we were designing a fleet of 
planes for an inexperienced company in a brand new market. Rather 
than spread the Prime Mover fleet thin throughout all of Aeroworld, 

the fleet was designed to serve a restricted area, the Northern 

Hemisphere. By doing this, a more reliable and efficient service could 
be provided. Concentrating on a restricted area and therefore 

providing an extremely dependable service would install a “piece of 
mind” among our customers. This eventually would attract more 
customers. From a marketing and economic point of view, this was 
deemed to be more important than serving a larger area at first and 
decreasing the efficiency of our service. The possibility of expanding 
the service throughout the whole world still exists, but as a new 
company creating a good reputation is of upmost importance. With this 
strategy in mind, the economic analysis could begin. 

One the most important results obtained from the economic 

analysis was that fuel costs, using an average price rate of $12.50 per 
milli-amp hour, accounted for 73% of the total fleet life cost. This 
differs quite a bit from the original design objectives. This was simply 
due to the fact the fuel cost equation was not interpreted correctly. 
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Figure 12.1 



0 Fuel Costs - 73 % 

M Production Costs - 25% 
M Maintenace Costs - 1 .5% 
H Operation Costs - 0.5% 


This result had a huge impact on the mission and final design. 
This large percentage drove Group Exodus to design optimally for a fleet 
that would serve a restricted area, namely the Northern Hemisphere of 
Aeroworld. The fuel costs also had repercussions on the final design. 
As can be seen in Figure 12.2, a one amp increase in the current draw at 
cruise results in an approximate 15% increase in the fleet life cost per 
volume moved (FCPVOL). 

Figure 12.2 


FLEET LIFE COST PER VOLUME MOVED 
vs . CURRENT DRAW 



FCPVOL 
($/cu. In.) 


76 


This translates into a more expensive price for our customers. 

Consequently, it was necessary to keep the current draw as low as 
possible. It was observed that one of the best ways to accomplish this 
was to minimize the weight of the aircraft. Therefore, Group Exodus 
opted for a relatively small cargo volume compared to other groups. 

The one drawback of this size volume was that the fleet size had to be 
increased; but, as mentioned before, this larger fleet size would insure 
the reliability and efficiency that it necessary for good overnight 
package delivery service. 

The pricing scheme adopted for this overnight package delivery 
service is quite simple. Customers will be charged per volume of the 
parcel with and added overseas charge for intercontinental shipping. 
There are two reasons for this overseas charge. The main reason is that 
after analyzing the daily cargo load for the Northern Hemisphere 

(Appendix I), it was noticed that 58% of this cargo is shipped overseas. 
Another reason for this overseas charge, is that customers will expect to 
pay more for parcels that are shipped overseas. Consequently, overseas 
shipping proved to be the best avenue for profit. In order to break 

even on the original investment in half of the fleet’s life, 150 days, a 
price of $8.74 per cubic inch for intracontinental shipping and a price 
of $11.01 per cubic inch for overseas shipping will be charged. The 
total fleet life cost is $57,069,600. Other pertinent results of the 
economic analysis are listed on the spread sheet in Appendix I. 

12.1 Production Costs 

As can be seen in the cost breakdown, the production cost 
accounts for 25% of the total fleet life cost, which makes it the second 
most important cost. Consequently, keeping this cost as low as possible 
translates into sizable price reductions for our customers. The 
production cost influenced the overall design of the Prime Mover. After 
reviewing the equation for the production cost, it was realized that the 
number of man hours was the most controllable aspects of the total 
production cost. Consequently, Group Exodus chose to keep the overall 
design as simple as possible in order, to keep the cost as low as possible. 
After reading past reports, it was estimated that 120 man-hours would 
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be sufficient to build the Prime Mover prototype, 
were estimated as follows: 


Construction costs 


CONSTRUCTION COSTS 


PROPULSION COST: 

Propeller (Zinger 11-5) 

$3.00 


Engine( Astro 15) 

$125.00 


Speed Controller 

$70.00 


Batteriesn21 

$15.00 



$213.00 

CONTROLS COST: Transmitter 



Receiver 
Battery Pack 
2 Servos 




$300.00 

STRUCTURE COST: 

Wood (balsa and spruce) 

$45.00 


Monokote 

$30.00 


Hinges, Clamps, Screws 

$10.00 


Landing Gear 

$22.00 


Glue. Tape 

$20.00 


$127.00 


TOTAL $640.00 

Using the given equation for production costs, a value of $376,000 
was obtained. 

12.2 Maintenance Costs 

The cost analysis (Figure 12.1) revealed that the maintenance 
costs only accounted for 1.5% of the total fleet life cost. As a result, this 
cost will not have a big impact on the customer. Nevertheless, it was 
minimized as much as possible. Maintenance costs are based on the 
time it takes to replace the batteries. Group Exodus therefore designed 
a removeable portion of the fuselage so that the batteries could be 
replaced as quickly as possible. It was estimated that a battery change 
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will take one minute and that the batteries will have to be changed on 
the average every 2 flights. Therefore, the maintenance costs were 
calculated to be $25 per flight. 

12.3 Operation and Fuel Costs 

Operation costs were estimated to be only 0.5 % of the total fleet 
life cost (Figure 12.1). Along with the maintenance costs, the operation 
costs have little impact on the customers. The operation costs are based 
the flight time in minutes at the maximum range and the number of 
servos in the aircraft. As a result of the simple design of the Prime 
Mover , the number of servos was limited to 2. The maximum flight time 
in minutes for our designed mission with a one minute loiter included is 
4.5 minutes. Consequently, the operation costs total a mere $9.00 per 
flight. 

Fuel costs as stated before, accounted for 73% of the total fleet life 
cost and therefore had the biggest impact on the customers. As a result 
these costs were kept to minimum as much as possible. This was 
described in the introduction of this section. Using an average value of 
$12.50 per milli-amp hour, the fuel costs were calculated to be $1,604 

per flight. 
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13.0 Results of Technology Demonstrator 

Development 

13.1 Configurational Data, Geometry, Weights and C.G. 

The configurational and geometry of the Prime Mover remained 
unchanged from the initial design concept. However, the methods of 
connecting the wing, tail and landing gear were altered. In addition, 
weight and C.G. location varied with respect to the predicted values. 

The wing was initially designed to be connected by hooks to a 
piece of spruce glued on top of the fuselage. However, Group Exodus 
desired to connect the wing in a more secure manner. A 3/16 inch 
piece of plywood was used to form a bulkhead inside the aircraft with 
an additional 1.5 inches on top for two holes. The wing had two birch 
dowel pegs that fit into these holes. The trailing edge was connected in 
the original manner. In future versions, the fuselage should be made 

wider so that a 4X4 inch square may be cut out to provide a continuous, 
uninterrupted cargo space. The tail was originally to be attached by a 
screw and hook design. The hook was replaced by velcro because it is a 

simpler technique. 

The main landing gear was placed further aft than originally 
planned. This was done to improve ground handling under the 

advisement of an expert in the field. Also, the main gear was bolted to 
the fuselage instead of being glued. This was proven to be an important 
decision, as several other aircraft with landing gear that did not use 
bolts were torn off. In future versions, the fuselage should be made 

slightly higher to account for thickness of the plywood support and the 
length of the bolt on the inside of the aircraft (approximately 1/4 inch). 

The avionics compartment should be changed in future versions. 
Two inches of available length in the fuselage may be used for the 

structure of a sturdier compartment. The present compartment is 1/16 
inch balsa and is very fragile. It is also a very tight fit and adds 

difficulty to installation. 

The batteries Group Exodus desired were not available, so lighter 
batteries with less power were used. Our motor was slightly lighter 

than noted on the specification sheet. The reason for the deviation is 

unknown. The main landing gear was slightly heavier than expected 
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because it was purchased instead of built. This was done to increase 
survivability and was probably one of the most important decisions 
made since the Prime Mover had several hard landings. The fuselage 
was heavier than expected due to the change in the method used for 
connecting the wing. Also, battery wire, control cables, and the avionics 
compartment were not accounted for in the original estimation. 

The C.G. was off by 0.5 inches. This was due to a combination of 
several small weight differences and the relocation of the main landing 
gear. To correct this, a 3 oz weight was placed in the rear of the 
aircraft. 

13.2 Flight Test Plan and Test Safety Considerations 

The flight test plan for the Prime Mover prototype consisted of an 
indoor flight in Loftus Center. The plane was designed to take-off 
within 20 yards and to safely complete a figure eight flight 
configuration within an area with a 100 yard length and a 40 yard 
width. In order to ensure a safe flight. Group Exodus meticulously 
completed a pre-flight check. Some of the major checks included: 
shaking the prototype to check for loose parts and the overall integrity 
of the aircraft, dropping the prototype from a height of three feet to 
check the soundness of the structure especially the landing gear and 

polyhedral wing joints, and testing the radio control and the 
corresponding movement of the control surfaces. 

13.3 Flight Test Results-Taxi and Controlled Flight Tests 

The taxi test was, for the most part, a success since it was 

confirmed that the Prime Mover prototype does indeed fly and that 

there were no major mechanical problems. Yet, the design group was 

alerted to a few minor problems that had to be attended to before the 
flight test. As the Prime Mover began its take-off, a bias in the rear 
landing gear caused the aircraft to veer to the left. The design group 
attributed this problem to the asymmetry of the rear wheel axle, which 
is curved on one side to clear the wheel. For flight tests, the gear was 
tuned to eliminate the bias and the axles were greased to decrease the 
force on the axle, but with little success. In retrospect, since the freely 
rotating rear axle does not provide adequate ground handling, it is 
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apparent that the rear gear should be attached to the rudder server 

Another factor that contributed to the tendency of t e pro o yp 
left during the taxi test was the warped rudder. Th,s warp resulted 

from the monokote clinging too tightly to one side, and therefo e 

creating a bend in the main spar. In an attempt.ng to correct th.s 
problem, the heating iron was applied to one side only to try to remo 
L warp. This attempt was unsuccessful. A better soluuon would be to 
make Z vertical stabilizer with suffer material. This would minimize 
the need for precise monoko.ing which is difficult for amateur but ders^ 
Finally, the taxi test revealed an inadequacy in the wing l 01 " ts - 
brief aborted flight, the aircraft made a hard landing, which cracked 

the spar webs that hold the wing joints m place. These we s were no 

designed to take such an impactive force, which occurred because a gap 
between the leading edges of the wing sections allowed a small amoun 
of wing rotation in the yaw axis. To strengthen these joints, the ongina , 
1/16 inch thick webs were replaced with 1/4 inch thick webs with the 

grain oriented vertically to prevent shear. During fhght lesls ’ t e 
between the wing sections were filled with small, pliable rubber shock 
absorbers to prevent transmission of impactive forces into t e joints. 

During the controlled flight test, a few additional problems were 
revealed in the construction of the Prime Mover. In flight, the aircra t 
tended to roll to the left, reducing controllability during right turns. 
The pilot believed that this was this result of slight wing warp, causing 
a roll moment during flight. This is the most likely cause of this 
problem, since the aircraft was balanced along the roll axis. no er 
possibility is that the rudder was in the wake of the outboard wing 
sections which are at 17“ actual dihedral. This was not accounted for in 
our calculations. Future analysis may show the need for a larger rudder 
,o compensate for this interference. The inability of the freely rotating 
rear gear to provide directional stability once again made ground 
handling difficult. Finally, another hard landing, this time on the nose 
of the aircraft, cracked the 1/4 inch thick balsa engine mount. Future 
versions should utilize a more durable material, such as spruce or birch 
plywood, for this critical structural member. 



13.4 Manufacturing and Cost Details 

The Exodus DR&O set a limit price of two hundred dollars for 
construction materials. Exodus believed that the technology 
demonstrator could be constructed for approximately one hundred and 
twenty five dollars. This is the figure used in obtaining the initial cargo 
fee. Though cost overrun kept Exodus from meeting the one hundred 
and twenty five dollar goal, the DR&O goal was accomplished. Final 
expenses totaled to $164.22. Figure 13.1 shows the percentage costs 
breakdown. A more detailed cost per component break down in shown 
in Figure 13.3. 

Major sources of pricing error included the monokote, landing 
gear, and excess wood. Monokoting the wing required two rolls of 
monokote instead of one. This increased the monokote price by $10.00. 
Initial estimates on the price of light weight wheels were off by 30%. 
Extra expenses in construction materials (balsa and spruce) arose due to 
manufacturing errors. In continued construction of the Exodus fleet, 
these additional costs can easily be avoided. 

The number of construction manhours spent on the Prime Mover 
prototype was 112. This is within 7% of the initial 120 manhour 
estimate. Again, this value should be reduced as manufacturing 
experience is gained. Figure 13.2 shows the component division of 
manhours 
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Figure 13.1 

Cost Breakdown by Aircraft Component 



■ fuselage 
H vertical tail 

■ horizontal tail 
[3 wing 

□ landing gear 

■ monokote 
EH control unit 
El propulsion 

□ misc 

□ tax 


Figure 13.2 

Manhours per Component 



Monokoting 

Empennage 

Fuselage 

Wing 

Controls 

Landing Ge 

Instrument I 

Propulsion 

Misceilanec 


15 . 89 % 
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IIPiiillilliiJM* , , w . . . . 

l/8"xl /8"x36" spruce diagonal beams (sides) 




1/8"x3"x36" balsa 
1 /4"x3"x36" balsa 
1 /4"x1 /4"x36" balsa 
1/8"x1 /4"x36" balsa 
1/1 6"x1 "x36" balsa 
Fuselage Total 

^'v '■■■■ !■.. ■ > - ‘ ' s - > , -.% J :• v. x <<• 

i/8 x 1/4 x 36" balsa 
1/8 x 1/2 x 36" balsa 
1/4" x 1/4" x36" balsa 
1/8 x 3 x 36 balsa 
Vertical Tail Total 


instrument prot wall 

firewall 

main frame 

vertical beams and misc 
cargo support panel 

15.51 


1.73 

2.4 

3.36 

4 

0.42 


balsa Internal 

balsa/tria Trailing Edge 
balsa Frame 

balsa block Attaching block 

3.38 

, -< - .;-Wc.sL 


wing mounting kit 

1/4"x1/4"x36" balsa Frame 

1/8"x1/4"x36" balsa Internal 

1/8 x 1/2 x 36" balsa/tria Trailing Edge 
Ho rizo n ta I Tail Total . 


CvfBsW; 

0.4 

0.69 

0.56 

1.73 


2.29 

0.56 

0.8 

0.69 


4.34 


1 

Birch Dow 


Wing connection 

1 

24"x6"x1/4" 

plywood 

polyhedral joints 

1 

1 2"x6"x3/1 6" 

plywood 

mounting bulkhead 

2 

1 / 4"x1 /4"x48" 

spruce 

main wing spars 

1 

1 / 2 ” x1/8"x48" 

spruce 

main wing T.E. 

3 

3/16x3/16x48" 

spruce 

all wing L.E. 

4 

1/4“x1/4"x36“ 

spruce 

side wing spars 

2 

1/2"x1/8"x36" 

balsa 

side wing T.E. 

8 

3"x36"x1/16” 

balsa 

ribs 

1 

1 /8"x3"x36" 

balsa 

ribs 


Wing Total 




25.18 


Landing Gear, 

1 Main Gear 

1 Tail Gear 

2 wheels 

1 tail wheel 

2 Collars 

Landing Gear Total 
Monokote |||i| lliS®. 

3 Clear monokote 

1 trans. blue monokote 
Monokote Total^^^ 

Control 

1 48" Push Rod 
1 36" Push Rod 
1 Control Horn Pack 
1 nylon hinge tabs 
Control Unit Total 
Propulsion::|f|^^ 

1 12x4 prop 

1 11x5 prop 

2 Velcro mounting 
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Appendix A 

Selection of Aspect Ratio 


The drag polar was determined by calculating the drag 
coefficient for various lift coefficients. The drag was divided into 

parasite drag and induced drag in the equation C D = C D o + Cl h teAK. 
The drag polars for the various aspect ratios are illustrated in Figure 
A 1 For a given lift coefficient, the drag coefficient decreases with 
increasing aspect ratio. Thus, it is desirable for the wing to have a 
large aspect ratio so that the drag will be minimal. A larger aspect 
ratio decreases the induced drag, the greatest contributor to the total 
drag for large lift coefficients. 


Figure A.l 

Drag Polars for Various Aspect Ratios 



The data from the drag polar was then used to generate the lift- 
to-drag ratios for the various lift coefficients as shown in Figure A.2. 
The exact value of the maximum lift-to-drag ratio was determined 
for various aspect ratios using the relation L/D) ma x = 
(C Do xeAR) 1/2 /2CDo and the corresponding lift coefficient was found 
using C L [atL/D)max] = SQRT(rce ARCdo) (Ref 1). As the aspect ratio 
varies, the angle of attack for zero lift remains constant at -2 degrees 
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and the stall angle occurs at an estimated 14 degrees. Tte » e« cency 
factor was assumed to vary with aspect ratio according to the 
equation e=l. 78 ( 1 . 00-0.045 ARO.«)-0.64 (Ref 3) -to- 

drag ratio increases with increasing aspect ratio. The * 
approximately 20% between an aspect ratio of 7 and LL The 

corresponding lift coefficients vary by approximately the ame 

, f ___. trc tn 7 At cruise conditions. Cl must oe 

nprcent and ranee from .oj to ./. m wuw '' 

greater than or equal to .64. It is desired that the auplane cruise MB 
L/D)max for greatest efficiency, however. L/D var,« only by^few 
nercent near the maximum value. For example, /D y 
approximately 1% for C L values between .5 and .7 for an , asp«t rmo 
Of 9. Thus, a large aspect ratio will produce the greatest 

aerodynamic results. 


Figure A.2 

Variation of Lift-to-Drag Ratio with Lift Coefficient 

w- Aenart R fltlOS 



n AR-7 

• AR-9 
■ AR-11 

• AR-13 


CL 


The effect of varying the aspect ratio on the lift curve slope is 
illustrated in Figure A.3. As the aspect ratio increases, the corrected 
lift curve slope increases. Thus for a given Cl the angle of Jt 
decreases with increasing aspect ratio whtch - — • 
the decrease is only minimal. For CLcruise equal to . , * 

attack will fall between 6 and 7 degrees for the range of aspect 
ratios. Therefore, aspect ratio does not have a great e 

reducing the required incidence angle. 


Figure A.3 

Lift Curve Slopes for Various Aspect Ratios 
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In order to assess the relation between the weight of the wing 
structure and the aspect ratio, the wing was simplified to a cantilever 
beam representing the main spar. This is a fairly accurate method 

for estimating the reaction of the wing to lift forces since the main 
spar is the primary carrier of stress due to bending about the roll 
axis. However, the spar weight for each aspect ratio is on y a 
characteristic weight for the wing. Assuming that the excess weig 
in the relatively ineffective middle of the spar is approximately 



equal to the weight of the balsa ribs and leading and trailing edge 
spars, the total wing weight is approximately equal to the spar 

weight plus about 5 ounces of monokote. 

For each aspect ratio, AR=b 2 /S, since S=bc is fixed, the span, b, 
and the chord length, c, are defined. In addition, since the spar 
height, h, is .117c, this is also determined by the AR. The spar is 
located at the thickest part of the airfoil since the maximum stress in 
a symmetric beam, a max , equals Mh/2I, where M is the bending 
moment and I is the moment of inertia (Ref. 8, p.116). Since the 
wing from root to tip is effectively a cantilever beam under an even 
force distribution, the moment equals Lb/8, where L is the total lift 
force (Ref. 8, p.162). For ease of calculation, the beam cross-section 
was modelled as a rectangle of height, h, and thickness, t, making 
I=th 3 /12 (Ref. 8, p.350). For a given aspect ratio and material with 
known yield stress, a spar thickness is defined. Given a matenal 
density, p, the weight of the spar is ptbh. 

Using an integral beam-bending analysis, the deflection of 
the wingtip is found to be Lb 3 /16EI, where E is the Youngs 
Modulus (Ref. 8, p.163). For spruce, E=1.3e6 psi, p=.01 61b/in 3 , 

and Otxiax = ^^^® psi. 


Figure A.4 


Spruce Wing Analysis Based on Main Spar Definition 
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Appendix B 

Calculation of the Efficiency Factor 


The efficiency factor was estimated by taking the average of 
the efficiency factor values obtained from the Empirical Formula and 

Jensen’s methods. 

1. Empirical Formula e=1.78*(l-0.045AR U(i8 )-0.64 

This equation, based on empirical data for straight-winged 
airplane, provides a good approximation although it was not 
determined for RPVs specifically (Ref 3). The efficiency factor is .765 
for an aspect ratio of 9.62. 


2. Jensen’s Method 1/e = l/ewing + 1/ebody + l/e 0 .l.cr 

In this equation, e w i n g and et,otly varying with aspect ratio an 
e„,hcr is approximated as 20 (Ref 3). Referring to Figure B1 .he 
value of ewing for a rectangular wing with an aspect ratio o ,s 

approximately .76. 

Figure B.l 

Wing Efficiency Factor 



To determine a value for e bo dy. a body efficiency parameter is 
defined as Ebody = CbodySbody/Srcf w * iere ^ body ' s l * ie cross 




area of the body and Sref is the wing area. The value of 1.22 for 
Ebody is obtained from Figure B.2 for a rectangular fuselage and an 
aspect ratio of 9.62. For S bo dy=20.23 inland S rc f=MI6 in . e bo d y is 
85.39. Thus, the value of the efficiency factor obtained from Jensens 

method is .726. 

Figure B.2 

Fuselage Efficiency Parameter 



The two values differ by only 5.1%. The average of these 
efficiency factor values is .746. This averaged value is used in all 
calculations containing the efficiency factor. 




Appendix C 
Drag Breakdown 


C 


The drag of the Prime Mover was predicted using the drag 
breakdown method outlined in Jensen’s thesis (Ref 3). This method 
estimates the total parasite drag by summing the contribution of each 
individual component of the airplane based on the wetted area and 

referencing it to the area of the wing according to the equation: 

c* - y r c /s r**f 


where Cf=.0055 was selected based on skin friction coefficients of other 
model aircraft. 


1. Fuselage 

S wet-945 in 2 
C Do =-003671 


2. Battery Pack 

S W et=35.6 in 2 
Cdo=-000138 


3. Frontal Section 
S we t=55.9 in 2 
C D o=-000217 


4. Rear Section 

S W et= 33.4 in 2 
C Do =000130 









5 . 


6 . 


Wing 

S we t= 2236 in 2 
Cdo=-008685 



S we t= 416 in 2 
Cdo=.001616 




7. 


8 . 


Vertical Stabilizer 
S W ei= 198 in 2 
C D o=-000769 



Landing Gear 
S W et=79.9 in 2 




The total wetted area for the aircraft is 3999.8 in 2 and, given a 
reference area of 1416 in 2 , the total parasite drag is .0155. Adding an 
additional 15% for interference, roughness, and proturberances provides a 
final estimated value for the parasite drag of .0179. 






Appendix D 
Sample Propeller Data 


A) 

B) 

C) 



Propeller Designation. 

Number of Blades: 2 

Select one of the following airfoil 


( Inches) 


D) 

E) 

F) 

G) 

H) 
I ) 


I-- ROP 

11-5 

Di ameter : 
sections: 

I NV I SC ID FLAT PLATE 
THIN FLAT PLATE 
SYMMETRICAL 
CLARK Y 
RAF -6 

NACA44X XLOWRE 
Blade thickness may be entered as either. 

1 ) Fraction of chord 
— v 2) Inches 

Bla de data may be entered at radiai 1 —ticn- specked ee 
1) Fractional Radius 
— > 2) Inches 

Radius* at wh i ch blade setting is measured: 

Blade setting (i.e. ref angle for mhole blade): 


1) 

2 ) 

3) 

4) 

5) 

6 ) 


(3-9) : 


14.86 

9 


Chord 

.707 
. 7B 
.873 
.924 
. 932 
.911 
. 843 
.756 
.632 


Thickness* 

Angl e 

.234 

38.51 

. 198 

27.95 

. 191 

21 . 7 

. 181 

17.66 

. 171 

14.86 

. 157 

12.81 

. 138 

11.25 

. 116 

10.02 

.095 

9.04 


J> 


K) 


L) 

M) 

N) 


Data Point Radius* 

l: 1 
2 : 1 .5 
3: 2 
4:2.5 
5:3 
6:3.5 
7:4 
8 : 4.5 
9:5 

Select desired refinment of analysis. + 

1) Analysis by simple blade element theory. 

2) Analysis including induced velocity. losses . 

> 3 ) Analysis including induced veloci y 

These Cl/Cd coefficient adjustments may be selected 

1) No Cl/Cd adjustments 

2) Mach number adjustment 

3) Reynolds number adjustment 

y 4 ) Mach and reynolds number adjustmen s 

Select altitude in thousands of feet. .0*-u 

Specify one of the foil owi ng : MPH 

> i) Airspeed FIXED at. 

2) Propeller RPM FIXED at: 

Range of Advance Ratio to be used in calculations: 

J min: . 2 J max : . 7 

NOTES: First data point nust be less than 3 il radius; others .ust progress outward. 
Designation oust start with a letter and »ay not contain a com.. 

For square tip blades (only) use tip as last data point. 

Angles lust be specified in degrees, lengths in inches, 
t Units s ust be as specified in lines D and E. 




cn cn 


Angles of Attack (Degrees) 


.§73 .066 .06 .053 .046 .038 .029 .021 .012 3E-03 -6E-03 

.028 .028 .027 .026 .024 .021 .018 .015 •«» ,E ' , _ 3 

.519 .582 .64 .683 .712 .737 .727 .6B9 .6 3 .264 -3.999 

.244 .226 .208 .194 .18 .168 .158 .149 .141 .33 .127 

9120 7514 6389 5557 4917 4409 3996 3654 3366 3120 2907 
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Appendix E 



Sample 

Cruise Conditions 

; for Electric Motor Performance 

input 

Name 

output Unit 

Comment 

Q 

1.0701 psf 

dynamic pressur 

.002378 

rho 

slug/ft3 

air density 

30 

vel 

ft/sec 

air speed 

Cd 

.03641082 

a/c drag coefficient 

. 0179 

Cdo 

zero lift drag coefficient 


Cl 

.63938933 

a/c lift coefficient 

• 74 

eff 


efficiency factor 

9.5 

AR 


aspect ratio 

1 

n 


load factor 

6.5 

W 

lb 

a/c weight 

Q 5 

s 

ft-ft 

wing area 


Preq 

15.05493 W 

a/c power required - level flight 


ROC 

-.0004831 ft/s 

rate of climb 


Pavail 

15.050672 W 

power available from propeller 


V 

9.2099145 volt 

armature voltage 

9.563 

vset 

volt 

battery voltage 

.08 

Kb 


battery constant 


i 

4.4135691 amp 

motor current draw 


motrpm 

11128.572 rpm 

motor speed (rpm) 

.12 

Ra 

ohm 

armature resistance 

.00078 

Kv 

volt/rpm 

motor speed constant 


proprps 

4675.8706 rpm 

propeller speed (rps) 

2.38 

gr 


gear ratio 


J 

.41979831 

propeller advance ratio 

.917 

propd 

ft 

propeller diameter 

eta 

.74524122 

propeller efficiency 

1.08 

Kt 

in-oz/amp 

motor torque constant 

.95 

greff 


gear efficiency 


fltime 

815.66639 sec 

flight time 

1 

batcap 

amp- hr 

battery capacity 


range 

CT 

CP 

24469.992 ft 

.0362363 

.0204121 

range 


moteff 

.54803682 

motor efficiency 

Rule 

oeff 

.40841963 

overall efficiency 


Q= . 5*rho*vel~2 
Cd=Cdo+C1^2/ (PI () *eff*AR) 

Cl= (n*W) / (Q*S) 

Preq=Q*S*Cd*vel 
ROC= (Pavail-Preq) /W 
v=vset-Kb*i 
mot rpm = (v-i*Ra) /Kv 
proprps=motrpm/ (60*gr) 

J=vel/ (proprps*propd) 

CT=Ct ( J) 

CP=Cp ( J) 

eta=Ct ( J) *J/Cp(J) 

Pavail=eta*Cp ( J) *rho*proprps A 3*propd A 5 

Cp ( J) *rho*proprps A 3*propd A 5= ( (Kt/Kv) * (v*i-i*i*Ra) * . 0005454-floss (motrpm) ) *gre 

f ltime=batcap/i 

range=vel*fltime*3600 

k=-9. 4225+6. 0529E-3*60*proprps-l . 3867E-6* (60*proprps) ~2 
f= (60*proprps) ~3*1 . 4337E-10- (60*proprps) A 4*5 . 5867E-15 
mat«ff“k+f 


Appendix F 

Equations and Assumptions 


F'l 


F.l Wing contribution to Cm 

Cm ow = Cm acw + Cl ow (Xcg/c - Xac/c) 

Cm aw = Cl a w (Xcg- Xac) 

assumptions: small angles - 15° max 

cos 15° = 1 (actually = .966) 

sin 15° = .262 (actually = .259) 

L/D » 1 (actually = 14) 

Zcg is negligable since multiplied by sine 


F.2 Horizontal stabilizer contribution to Cm 
Cm 0 h = ^h* v h* cl och*( £ o + J w - h ) 

Cm ah = -n h *V h *Cl ah *(l - d£/ da) 

assumptions: = $1 this is an estimate using ref.6 pg 47 

£ = £ 0 + (d£/ da) a w } 

£ 0 = 2*Cl ow /(ftAR w ) } ref.6 pg 47-48 

d£/ da = 2*Cl aw /(7lAR w ) ) 
lt/Zcgt » 1 (actually = 24) 


Fuselage 

contribution 

to 

Cm 

X 

< 

+ 

Cm 0 f = 

k2 -kl 

_Z 

wf 2 (a ow 


36.5 *S w *c w 




Cm a f = 

1 

Z 

wf 2 (d£/ da) Ax 


36.5 *S w *c w 




assumptions: k2-kl 

= 

0.95 


station 

Ax 


wf 

aow+if 

1 

0.33 


0.283 

-16 

2 

1.83 


0.396 

-2 

3 

1.48 


0.396 

-2 

4 

0.27 


0.313 

+32 


ref.6 pg 49-51 


wf 2 (aow+if)Ax 

-.427 
-.575 
-.464 
+.849 
total = -.616 



station 

Ax 

wf 

1 

0.33 

0.283 

2 

1.83 

0.396 

3 

1.48 

0.396 

4 

0.27 

0.313 


X 

d £/ da 

w f 2 ( d8/ da )Ax 

1.99 

1.10 

0.03 

0.92 

1.20 

0.344 

0.74 

0.093 

0.022 

1.62 

0.08 

0.002 
total = 0.398 


F.4 Setting the wing position 

Meg = Mac + Lw (Xcg - Xac) 

Cmcg = Cmac w + Cl w (Xcg/c - Xac/c) 

Xac/c = Xcg/c + (Cmac w /Cl w ) 

F.5 Wing and Fuselage contribution to Cn 
Cn Bw f = *K n Krl — fs— f 
S w ^ 

assumptions: K n = .00075 

Kri = 1-4 

F.6 Vertical stabilizer contirbution to Cn 
Cn Bv = V v Cl av Tlv(^ d °/ d P) 

assumptions: 1 V (1+ da/ dP) = 0.992 + 1.53 (S v /S w ) + 0 009 AR w 

ref. 6 pg 71 


F.7 C1 B calculations (ref.4) 

(Cl B ) t0 tal = ( C1 b)w + (Cl B )v + (AC1 b )i + (AC1 b ) 2 + (C1 B )w. r=0 
assumptions: 

(Cl B )w = - 20799 

(CI B )v =-.01129 

(AC1 b )i = -.0006 

(AC1 b ) 2 =+.00016 

(CI B )w. r =0 = --00364 


F.8 Required EDA and actual panel angle(ref.4) 


30 ft/s => 22°/s roll rate 
desire 25°/s roll rate 
£ =25 * 1 16.16 

10 ° 22 100 
T = 13.2° 


Roll Rote 
(deg/sec) 



10 20 30 40 50 60 70 60 90 100 

Airspeed 
(ft/sec) 


Break Point = 40% span 
Moment Fraction = 0.775 
13.2° EDA = 0.775 * panel angle 
Panel angle = 17° 

w Fraction 



p 9 Required Yaw Angle (ref.4) 


o o i 
root 


0 2 0 3 0 4 0 5 0 6 0 7 

Semi-span Station 


00 


0.9 I 
tlO 


Roll Rate 
(deg/ sec) 



25°/s = 100 * 13.2 

X 116.16 10 


X = 22°/s 


0 2 4 


6 a 10 12 14 16 10 20 

Yaw Angle (degrees! 


Appendix G 

Performance Calculations 


The three computer codes listed use the same procedure to find 
a motor speed for a particular flight velocity. In "takeoff.f", lines 
through 22 input the necessary variables. Note that t e user mus 
manually input the propeller coefficients as indicated in ltnes 67 
through 78. In addition, a motor loss has been included to me 65 
using a motor efficiency as a function of propeller RPM. Ftna y, 
gearing efficiency is found in line 81 as the value 0.95. 

Lines 56 through 89 of "takeoff.f is the iterative process that 
finds the converging motor speed for a particular flight velocity. 
Once the motor speed has converged within a difference of 20 (line 
89) the program utilizes the values found for motor RPM, coefftcien 

of thrust and power, and current. _ 

After the motor convergence is complete, various ot er 

calculations can begin. For each particular flight velocity and voltage 
setting, values for power, thrust, lift, and drag can be computed, 
thereby allowing an analysis of the kinematics of the aircraft 
"Takeoff.f" uses a time-step process in lines 106 throug 
calculating distance travelled and velocity for each increment in 
time. When the velocity of the aircraft exceeds the takeoff velocity 
the program terminates. "Perf.f" and "pavail.f" both calculate the 
power available, power required, and rate of climb at each esire 
flight velocity and voltage setting. "Perf.f" finds the voltage setttng 
and current require for cruise conditions (rate of dim 
approximately zero) and then calculates the range and endurance for 

each flight velocity. 


Landing Calculations 


Xqr =W -Un[l + B-VtoI 

8 2B A 

A=nW 

B=C[*5S 


VTD=V sla n- 22 .6 ft/ S 


W=6.7 lb 

|i=0.15 


S=9.68 

C D =0.01 


9+ LJ Ai 

tc(.74)(9.68) 


Q-2 


0.077 


Xgr= 43 .8 ft 



Gliding Calculations 


G- 3 



Tmin=3.33° 

Glide Dist)max=L/D) max x altitude =17.2 x 25 ft = 430 ft 



takeoff . f 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 c 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 85 

57 

58 

59 

60 15 

61 
62 

63 

64 


real j , kt , kv, mass, mrpm,mrpmn, mrps, mu, lift 

open ( 10 , file = f take' ) 
read (10,*) wem 
read ( 10 f *) cargo 
read (10, *) rho 
read (10,*) clmax 
read ( 10 , * ) clto 
read (10,*) cdo 
read(10,*) sref 
read ( 10 , * ) e 
read ( 10 , * ) ar 
readdO,*) dia 
read(10,*) kt 
read(10,*) kv 
readdO,*) rarm 
readdO,*) rbat 
readdO,*) batcap 
readdO,*) fusamp 
readdO,*) gearat 
read(10,*) mu 
close (10) 

dt-.Ol 

bvolts=l 4 . 4 
dia^dia/ 12 . 

pi=-4 . *atan (-1 . ) 

dia4=dia**4 

dia5=dia**5 

rtot=rarm+rbat 

open ( 100 , file=' takeout' ) 

do 40 c=0 . , 32 » , 4 . 
w=wem+c/16 . 
mass-w/ 32 . 17 4 
clto=clmax/l . 44 
cdto=cdo+clto**2/ (pi*e*ar ) 
ampmax=bvolts / (2.*rtot) 

pomax= ( (kt*bvolts**2) / (4.*rtot*kv) ) *2.*px 

pomxhp=pomax*l . 578e-7 

pomxwt=pomxhp/ 1 . 341e-3 

f acl=bvolts/ (2 . *rtot) 

fac2= (bvolts/rtot ) **2 

batlos=0 . 

time=0 . 

icount=0 

v=0 . 

s-0 . 

vto=sqrt (2 . *w/ (rho*sref *clmax) ) *1.2 
amps=0 . 

mrpm= (bvolts-ampmax*rtot) /kv 

continue 

icount=icount+l 

iter=0 
continue 
mrps=mrpm/ 60 . 
iter=iter+l 
prps**mrps / gearat 
prpm=prps*60 . 



f _o J 


takeoff . f 

65 

5 *prpm* *4 

66 

67 c 

68 c 

69 c 

70 c 

71 

72 

73 c 

74 c 

75 c 

76 c 

77 c 

78 c 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 c 


etam=-9 . 4225+6 . 052 9e-3*prpm-l . 38 67e-6*prpm**2+l . 4337e 10*prpm 


**3-5.5867e-l 


j-v/ (prps*dia) 

ct= 1407 9-8 . 2493e-2* j-8 . 7552e-2* j **2-7 . 4503e-2* j**3 
cp=6.4266e-4+. 31747* j-. 99812* j **2+1. 5975* j**3-l. 11 67 *j**4 

ct=9 7 987e-2-.11367*j-5. 0432e-2* j **2-6 . 2163e-2* j **3 

cp=2. 7985e-2-3. 7795e-3* j+9 . 3626e-2* j**2- . 38324* j**3+. 20948 *j**4 

12-4 . „ 

ct=7 . 6014e-2-. 10035* j-8. 643e-2* j**2- . 10066*7**3 
cp=2 . 0033e-2+l . 566e-3* j+1 .7352e-2* j**2-.18709* j**3 
12-6 

ct=. 12 5 97-. 167 33* j-1 . 4419e-2* j**2-4 . 9687e-2* j**3 
cp=3.399e-2+2 . 6607e-2* j-6 . 2701e-2* j**2-7 . 5545e-2* j**3 
pmot=cp*rho*prps**3*dia5*l . 152e4/ etam 
if (pmot .ge. pomax) goto 1001 

amps=facl- . 5*sqrt (fac2-4 . *pmot/ . 95*kv/ (rtot*kt*2 . *pi) ) 

mrprnn= (bvolt s— amps * rt ot ) /kv 

d=abs (mrpm-mrpmn) 

mrpm 3s mrpm+ (mrpmn-mrpm) * . 5 

if (iter .gt. 1000) then 

print * f "motor speed calculation did not converge 

goto 40 

endif 

if (d .gt. 20.) goto 15 
if (amps .gt. fusamp) then 
print *, "fuse current exceeded" 
stop 
endif 

prps=mrpm/ (60.*gearat) 

thrust=ct *rho^prps**2*dia4 

drag= . 5*rho*v**2*sref *cdto 

lift*. 5*rho*v**2*sref *clto 

frict=mu* (w-lift) 

if (thrust .It. frict) then 

print *, "friction exceeds static thrust" 

stop 

endif 

accel= (thrust-drag-f rict ) /mass 


delv=accel*dt 

ds=v*dt 

dbat=dt *amps/3600 . 

v=v+delv 

s=s+ds 

bat los— bat los+dbat 

t ime=t ime+dt 

if (s .gt. 200.) then 

print *, "distance greater than 200 ft" 

goto 40 

endif 

if (time .gt. 60.) then 

print *, "time greater than 60 sec" 

stop 

endif 

if (v .It. vto) goto 85 

wl=w/ sref 

write (100, *) s 

roc= (thrust-drag) /w*v 

gamma~asin( (thrust-drag) /w) 

print *, c, s 

print * , roc 





takeoff . f 

128 c 

129 c 

130 c 

131 c 

132 c 

133 c 

134 c 

135 c 

136 c 

137 c 

138 c 

139 40 

140 

141 

142 

143 

144 1001 

145 

146 

147 


print *, "V takeoff =" , vto 
print *, "Time for run (sec) =", time 
print *, "V at TO (ft/sec) =",v 
print *, "Distance (ft)=",s 

print *, "Battery Drain (amp-hrs) =",batlos 
print *, "Advance Ratio at TO=", j 
print *, "Thrust (lb) at TO=", thrust 
print "Lift (lb) at TO = "# lift 
print *, "Drag (lb) at TO=", drag 
print *, "Friction (lb) at TO=", frict 
print *, "Current Draw at TO (amps)='\ amps 
continue 

close (100) 
stop 

print *, "Pmot exceeded Pmax" 

print *, "Pmot=" , pmot, " Pmax=", pomax 

stop 

end 


perf . f 


real j , kt , kv, mass, mrpm, mrpmn, mrps 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 c 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 85 

46 

47 

48 

49 

50 

5*prpm**4 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 c 

61 
62 
63 


open (10, file =' stuff' 

read (10 , * ) 

w 

read ( 10 , * ) 

cargo 

read(10, *) 

rho 

read ( 10 , *) 

clmax 

readdO, *) 

cdo 

read (10, *) 

sref 

read(10,*) 

e 

read ( 10 , * ) 

ar 

read ( 10 , * ) 

dia 

read ( 10 , * ) 

kt 

readdO,*) 

kv 

readdO, *) 

rarm 

read ( 10 , * ) 

rbat 

read ( 10 , * ) 

batcap 

readdO, *) 

fusamp 

read ( 10 , *) 
close (10) 

gearat 


pi=-4 . *atan (-1 . ) 
dia4=dia* * 4 
dia5=dia* *5 


mass=w/32 . 17 4 
rtot=rarm+rbat 
w-w+cargo/16 . 
vel=30 . 

open (100, file~'perfout' ) 
do 35 vel— 0 . , 60 . , 1 . 
cl=w/ ( .5*rho*vel**2*sref ) 
cd=cdo+cl**2/ (pi*e*ar) 
drg= . 5*rho*vel**2*sref *cd 
preq=drg* vel* 1 . 356 
roco=10 . 

do 60 bvolts=6 . , 16 . , . 005 
ampmax=b volts/ (2.*rtot) 

pomax= ( (kt *bvolts**2 ) / <4.*rtot*kv) ) *2 . *pi 

pomxhp=pomax* 1 . 57 8e-7 

pomxwt =pomxhp/ 1 . 341e-3 

facl-bvolts/ (2 . *rtot) 

f ac2= (bvolts/rtot) **2 

mrpm= (bvolt s-ampmax*rtot ) /kv 

iter=0 

continue 

mrps=mrpm/ 60 . 

iter=iter+l 


prps=mrps/ gearat 
prpm=prps* 60 . 

etam=-9 . 4225+6 . 052 9e-3*prpm-l 


. 38 67e-6*prpm**2 + l . 4337e-10*prpm**3-5 . 58 67e-l 


i=vel/ (prps*dia) 

ct=9 7987e-2-. 11367* j-5 . 0432e-2 * j** 2-6 . 2163e-2*]**3 

cp=2 *. 7985e-2-3 .7795e-3* j+9 . 3626e-2* j**2- . 38324* j**3+. 20948*^**4 

pmot=cp*rho*prps**3*dia5*1.152e4/etam 

amps=facl- . 5*sqrt (fac2-4 . * (pmot/ . 95) *kv/ (rtot*kt*2 . *pi) ) 

mrpmn= (bvolts-amps*rtot) /kv 

d=abs (mrpm-mrpmn ) 

mrpm=mrpm+ (mrpmn— mrpm) * . 5 

if (iter .gt. 1000) then 

print *, "motor speed calculation did not converge" 

goto 60 

endif 

if (d . gt‘. 20.) goto 85 


perf . f 


64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 60 

85 

86 35 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 


if (amps ,gt. fusamp) then 
print *, "fuse current exceeded" 
goto 60 
endif 

prps=mrpm/ (60.*gearat) 
thrust=ct *rho*prps**2*dia4 
pavail=thrust *vel 
pavlhp=pavail/550 . 
pavlwt=pavlhp/ 1 . 34 le-3 
roc= (pavlwt-preq) /w 
if (abs(roc) .gt. .1) then 
goto 60 

else 

if (abs(roc) .It. roco) then 
endur=batcap/ amps 
range=vel*endur *3600 . 
battloss=amps/ vel/3600 . 
roco=abs (roc) 
endif 

endif 

continue 

write (100,*) vel, range, endur, battloss 
continue 

close (100) 
stop 

print *, "Pmot exceeded Pmax" 

print *, "Pmot=” , pmot, " Pmax=",pomax 

stop 

end 


1001 



pavail . f 


real j , kt , kv, mass,mrpm, mrpmn, mrps 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 60 

36 c 

37 

38 

39 

40 

41 

42 

43 

44 

45 85 

46 

47 

48 

49 

50 

5*prpm** 4 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


open ( 10 , file =' stuff ' ) 

read ( 10, *) w 

read (10 , * ) rho 

read(10, *) clmax 

read (10 , *) cdo 

read (10 , *) sref 

read ( 10 , * ) e 

read ( 10 , *) ar 

read (10, *) dia 

read ( 10 , * ) kt 

read(10,*) kv 

read(10,*) rarm 

read{10,*) rbat 

read (10,*) batcap 

read (10,*) fusamp 

read (10, *) gearat 

close (10 ) 

bvolts=14 . 4 

pi=-4 . *atan (-1 . ) 

dia4=dia**4 

dia5=dia**5 

mass=w/ 32 . 174 

rtot-rarm+rbat 


open (100, file-' pavout' ) 
open (101, file-' rocout' ) 
open (102, file-'prout' ) 
do 35 vel~0 . , 70 . , 5 . 

c l=w / ( . 5*rho*vel**2*sref ) 
cd=cdo+cl**2/ (pi*e*ar ) 
drg- . 5*rho*vel**2*sref *cd 
preq=drg*vel*l . 356 
continue 

do 60 bvolts=7 . , 14 . 4, . 1 
ampmax=bvolts/ (2.*rtot) 

pomax= ( (kt*bvolts**2) / (4.*rtot*kv) ) *2.*pi 

pomxhp=pomax*l . 578e-7 

pomxwt =pomxhp/ 1 . 341e-3 

f acl=bvolts/ (2 . *rtot ) 

fac2= (bvolts/rtot) **2 

mrpm= (bvolts—ampmax*rtot ) /kv 

iter=0 

continue 

mrps— mrpm/ 60 . 

iter=iter+l 


prps— mrps/ gearat 
prpm-prps *60. 

etam— 9 . 4225+6 . 0529e-3*prpm-l 


3867e-6*prpm**2+1.4337e-10*prpm**3-5.5867e-l 


c t =9 1 7987e-2- 1 11367* j-5 .0432e-2* j**2-6.2163e-2* j**3 

cp=2 1 7985e-2-3 .7795e-3* j+9 . 3626e-2*j**2-.38324*j**3+ .20948*3**4 

pmot-cp*rho*prps**3*dia5*l . 152e4/etam 

if (pmot .ge. pomax) goto 1001 _ ^ 

amps-f acl- . 5*sqrt (f ac2-4 . * (pmot/ . 95) *kv/ (rtot*kt 2. pi)) 

mrpmn— (bvo Its -amps* rtot) /kv 
d-abs (mrpm-mrpmn) 
mrpm— mrpm+ (mrpmn-mrpm) * . 5 
if (iter .gt. 1000) then 

print *, "motor speed calculation did not converge 

goto 60 
endif 


pavail . f 


64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 c60 

77 

78 35 

79 

80 
81 
82 

83 

84 

85 1001 
86 

87 

88 


if (d .gt. 20.) goto 85 
if (amps .gt. fusamp) then 
print *, " fuse current exceeded" 

goto 60 
endif 

thrust=ct*rho*prps**2*dia4 

pavail=thrust*vel 

pavlhp=pavail/ 550 . 

pavlwt=*pavlhp/ 1 . 341e-3 

roc= (pavlwt-preq) /w 

write (100, *) vel, bvolts, pavlwt 

write(101,*) vel, bvolts, roc 

continue 

write (102,*) vel,preq 
continue 

close (100) 
close (101) 
close (102) 
stop 

print *, "Pmot exceeded Pmax" 

print *, "Pmot=" , pmot, " Pmax=" ,pomax 

stop 

end 


Appendix H 

Two-dimensional Truss Analysis 

The Prime Mover will be subjected to a variety of loads during 
it’s life. These vary from lifting forces, to weight forces, to catapult 
forces. Each of these forces offers their individual challenges in 

fuselage design. 

It was Exodus's goal to produce an aircraft fuselage capable o 
withstand a load factor of 4.0 in landing situations, 2.5 in cruise 
conditions, and 2.0 during catapult launching. Due to the complexity 
of structural testing, three design options were studied. Each of 
these were analyzed separately. 

By the assumption that the forces on the aircraft are nearly 
symmetric, one side of the fuselage is analyzed in a two dimensional 
truss method. Fixing certain junctures of the truss design allowed for 
a static analysis. 

The forces were simulated by point loads at fuselage nodal 
locations in a two dimensional truss analysis, where each node is 
defined by a junction of two or members. Therefore, if a force was to 
be distributed between nodes 5 and 6, half the load would be placed 
at each location. Forces included the following depending on the load 
environment: landing gear weight, equipment weight, cargo 

weight, wing weight, battery weight, engine weight, empennage 
weight, fuselage weight, wing lift force, landing forces, engine forces, 

and catapult forces. 

These forces were multiplied by the desired load factors. The 
optimal designs of each initial design were compared for given nodal 
locations. The saw tooth truss design was chosen as the strongest 
because of it's ability to handle compressive and tensional loads. 

The program was then altered to allow for further optimization. 
The next two pages include the input data file used to analyze the 
final concept in a 4.0 g landing environment. The program that 
follows determines the minimal cross sectional area of any member 
dependent on material properties and the internal forces produced 
from the given loads. 
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FUSELAGE-2- 


36 

23 

8 

37 

22 

9 

38 

21 

10 

39 

20 

11 

40 

19 

12 

41 

18 

13 

42 

17 

14 

43 

29 

3 

44 

3 

27 

45 

27 

5 

46 

5 

25 

47 

25 

7 

48 

7 

23 

49 

23 

9 

50 

9 

21 

51 

21 

11 

52 

11 

19 

53 

19 

13 

54 

13 

17 

55 

17 

15 


.03125 

.03125 

.03125 

.03125 

.03125 

.03125 

.03125 

.0156 

.0156 

.0156 

.0156 

.0156 

.0156 

.0156 

.0156 

.0156 

.0156 

.0156 

.0156 

.0156 


65000 . 

65000 . 

65000 . 

65000 . 

65000 . 

65000 . 

65000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 

1300000 . 


o o o o o o o 


c 

c 

C 


Two Dimensional Truss Analysis and Optimization Pro 9 r ^ 

Internal Stress and Nodal displacement were compared to 
that of SPACETRUSS and found extremely accurate. 

The initial areas of the members are read from the input file 
entitled Alanding. These areas are varied based on decreasi g 
tSe area of the member with the highest load factor without 
exceeded structural limitation of a 4.0 g landing load. 

Dimension node (30) ,memb (60) , xpos ^0) , ypos <^0) , zpos ( ^ 12 ^ 125^ X ( 125) , ( CIV(12: 

££ 

INTEGER CR, start (60) , finish (60) 

OPEN(12, file=' Alanding' ) ,, J _, ., 

WRITE (6,*) 'Fuselage design A at ground conditions n= 

PI-4 . *atan (1 . ) 

Density* . 0058 

READ (12, ' (a) ' ) junk 

WRITE (6,*) . , . , . „ , 

WRITE (6,*) ' This is a 2 dimensional analysis. , 

write ( 6 *) ' All nodes are constrained m the y component. 

SSt) ' This is an analysis of the side of the fuselage only.' 

READ (12#*) nnodes, nmembs, junkl, junk2 

WRITE (6,*) 'There are nnodes nodes and', nmembs, members 
N=2*nnodes 

THIS DO LOOP ALLOWS FOR THE READING OF THE POSITIONS OF EACH 
JOINT TO BE USED IN CALCULATING LENGTH, THETA, AND THE 
SPRING CONSTANT VALUES 


For=0 
moment =0 
Kount=0 


20 


S^ADU^'t' n^(la) , junk!, Junk2, jnnk3, «pos (la, , ypos <I * ) 


If 

If 


( junkl 
( junk3 


•eq. 
■ eq. 


0 ) 

0 ) 


WRITE (6, *) 
WRITE (6, *! 


Constrained in x 
'Constrained in z 


direction 

direction 


at 

at 


la 

node' , la 


For=For+zfor (la) 

moment=moment+xpos (la) *zfor (la) 


Continue 


WRITE (6, *) 

WRITE (6,*) 'Sum of the Forces -', For 
WRITE (6,*) 'Sum of the Moments =' , moment 
WRITE (6, *) 


789 Kount=Kount+l 
cop=0 
TOTWE=0 
DO 10 1=1, N 
DO 10 J=1,N 
CIV(I, J)=0 
10 C(I,J)=0 

D— 0 


“ (Koint 1 '.^ 6 ^) READ (12,*) memb (Ib) , start (lb) , finish (lb) , Area (lb) ,E (lb) 
If (ib .ge. 44) DENSITY=.016 
JA=start (Ib) 

JB=f inish (Ib) , 

D=SQRT ( (zpos (JB) -zpos (JA) ) **2 + (xpos(JB) xpos(JA)) ) 

THETA=ATAN2 (zpos ( JB) -zpos ( JA) , xpos ( JB) -xpos ( JA) ) 

KSPRING=E (Ib) *Area (Ib) /D 
WETGHT (Ib) =D*Area (Ib) *DENSITY 


o o 


TMOIN= (Area (lb) **2) /12 

TOTWE=TOTWE+WEIGHT ( lb) 

Per (lb) = PI*PI*E(Ib) *TMOIN/ (D**2) 

*******G 0 i n g to the subroutine************!!! 
NJOINTS=nnodes 

CALL KMAT (C, THETA, KSPRING, JA, JB,NJOINTS) 
0 *******Made it back! ! ! ! 


LENGTH (JA, JB)=D 
ANGLE ( JA, JB) =THETA 
KFORS ( JA, JB) =KSPRING 
30 Continue 


***************** PARTITION MATRIX *************** 

THE CRITICAL ROW NUMBER WILL BE READ- MEANING THE ROW 


IN WHICH DISPLACEMENT: 


CR-7 

NEWA=CR-1 

NEWN=N-CR+1 


Do 42 lit=l,N 
Do 43 lat=l,N 
CIV (lit, lat) =C (lit, lat ) 
c WRITE ( 6 , *) lit, lat, C (lit, lat) 

43 Continue 
42 Continue 

Do 761 lit=l , N 
CIV (2, lit) =C (29, lit) 

CIV (3, lit )=C( 30, lit) 

CIV (4, lit )=C( 31, lit) 

CIV (5, lit ) =C ( 32 , lit ) 

CIV (6, lit ) =C ( 58 , lit ) 

Do 762 lat-7,33 

762 CIV (lat, lit ) =C (lat-5, lit ) 

Do 763 lat=34 , 58 

763 CIV(lat, lit) =C (lat-1, lit) 

761 continue 

Do 764 lit=l, N 
C (lit, 2) =CIV (lit, 29) 

C (lit, 3) =CIV (lit, 30) 

C (lit , 4 ) =CIV (lit , 31) 

C (lit, 5)=CIV (lit, 32) 

C (lit, 6) “CIV (lit, 58 ) 

Do 765 lat=7, 33 

765 C (lit, lat ) =CIV (lit, lat-5) 

Do 766 lat=34,58 

766 C(lit, lat)=CIV(lit, lat-1) 

764 continue 

C Input forces 

Do 767 1=1 , NEWN 

767 FII(I)“0 

do 768 1=2,28,2 

768 FII (I-l)-zfor (1/2) 
do 769 1=30 , NEWN+1 , 2 

769 FII (I-l)=zfor (1/2+2) 

DO 60 Id=l , NEWN 
DO 60 Ie=l , NEWN 
CIV (Id, Ie ) =C ( Id+NEWA, Ie+NEWA) 
c Write ( 6 , *) Id,Ie, ' ',C(Id,Ie) 

60 Continue 





C**************** Gaussian Section 
DO 310 ib=l,N 
X (ib) =0 

310 continue 

Do 330 K=l, NEWN 
CC= CIV (K, K) 

FII (K)= FII (K) /CC 

DO 320 I=K, NEWN 
320 CIV (K, I ) = CIV (K, I ) /CC 

DO 330 I=K+1 , NEWN 
CC= CIV ( I , K) 

FII ( I ) = FII ( I ) - CC*FII (K) 

DO 330 J=K, NEWN 
CIV (I, J)=CIV(I, J) -CIV (K, J) *CC 
330 continue 

DO 340 Ih=NEWN, 1, -1 
SUM=0 

DO 350 I j=NEWN, Ih, -1 
PROD=CIV ( Ih, I j ) *X ( I j ) 

350 SUM=SUM+PROD 

340 X(Ih)=FII (Ih) -SUM 

C*********************DEFINE results*************************** 


DO 80 IA=57, 33,-1 
80 X (IA) =X ( I A— 5 ) 

DO 90 IB=28 , 2 , -1 
90 X (IB) =X (IB-1) 

X(1)=0. 

X (29) =0 . 

X (30) =0 . 

X (31) =0 . 

X (32) =0 . 

X (58) =0 . 

DO 93 INA=l,nnodes 
IOP=2*INA-l 
IUP=2*INA 
U (INA) =X ( IOP) 

93 V ( INA) =X ( IUP ) 

C *********************Compare these displacements to SPACETRUSS***** 

c DO 95 IYP=l,nnodes 

C WRITE ( 6 , * ) 'U' , IYP,U(IYP) 

c 95 WRITE ( 6 , * ) ' V' , IYP, V ( IYP) 


C*******That completes the part on displacements now we must find the reaction fore* 

DO 222 I=l,nnodes 
FX (I) =xfor (I) 

222 FY (I) -zfor (I) 


C****** ****** **USING THESE FORCES AND DISPLACEMENTS 



< l- r 


c 

c 

102 

C 


c 


c 

105 


Deflection 

DO 102 1=1, nnodes 
DEF (I)=SQRT (U (I) **2+V (I) **2) 
WRITE ( 6 , * ) ' DEF {' ,1 , DEF ( I ) 

Continue 


Forces and Stresses 

tackle=0 

tbone=0 


DO 105 J=l, nmembs 
KA=start ( J) 

KB=f inish (J) 


tat=0 

FORCE (KA, KB) = (KFORS (KA, KB) * ( (U (KB) -U (KA) 
If (FORCE (KA, KB) .It. 0) then 
Forrat= (FORCE (KA, KB) /-Per (J) ) 


) *COS (ANGLE (KA, KB) ) + (V (KB) -V (KA) 


else 

Forrat=l 


STRESS (KA, KB) =FORCE (KA, KB) /AREA ( J) 

SHEAR (J) “STRESS (KA, KB) /2 

If (STRESS (KA, KB) .gt. 0) then 

If(J It 44) Strrat=STRESS (KA, KB) /400 

If(J ! ge . 44) Strrat=STRESS (KA, KB) /750 


else 

Strrat=l 


endif 


★ J 


I£(FORCE(KA,KB) .eq. 0) Forrat-l Exc( , eded in membe r' , J, ' Are, 

IWTtYrll :S: 1! SStJU:*!- —Stress Ratio Exceeded in member' , J, ' «. 

SHEAR (j,/2o5 he qt. 1) write <6,., '-Shear Ratio Exceeded in memher'.J,' 

?f1sHEAR(J)/750 .gt. 1) Write (6, *) ' —Shear Ratio Exceeded in member', J,' 
endif 


Testing Force and Stress Ratios 


If (ABS (Forrat-l) .gt. data) then 

tackle=Forrat 

tbone=Area ( J) 

data=Abs (Forrat-l) 

recall=J 

endif 

If (Abs (Strrat-1 ) .gt. data) then 

tackle-Strrat 

tbone=Area ( J) 

data=Abs (Strrat-1) 

recall” J 

endif 


Safety Factor 


CONTINUE 


' Member*' , recall 

' Ratio-' , tackle 

• Current Area=',tbone 

'Would you like to change area?' 
' If yes, input member' 

'If no, input 999' 


Write (6, *) 
Write (6, *) 
Write (6, *) 
Write (6, *) 
Write (6, *) 
Write (6, *) 


Read (5, *) recall 

If (recall .gt. 990) goto 212 

Write (6,*) 'Input area' 

Read (5, *) exit 


912 


913 


914 


915 


212 

737 

747 


If ((recall . ge . 1) .and. 
do 912 i-1,14 
Area (i) =exit 
write (6, *) 'Area 
else if ((recall .ge. 16) 
do 913 i-16,28 

Area (i) =exit 
write(6,*) 'Area 
else if ((recall .ge. 30) 
do 914 i=30 , 4 2 
Area (i) -exit 
write (6,*) 'Area 
else if ((recall .ge. 43) 
do 915 i=43, 54 
Area (i) -exit 
write (6,*) 'Area 

else 

Area (recall) -exit 
write (6,*) 'Area 


(recall . le . 14)) then 
1-14' , Area(l) 

.and. (recall .le. 28)) then 

16-28' /Area(l6) 

.and. (recall .le. 42)) then 

30-42' , Area(30) 

.and. (recall .le. 54)) then 

43-54' ,Area(44) 

recall,' ',exit 


data=0 
exit=0 
endif 
Goto 789 


WRITE (6, *) 
Write ( 6 , * ) 
FORMAT (2x, 
FORMAT (2x, 


'The total weight of the structure in onces-' 

' Member' , 8x, ' TENSION' , 9x, ' STRESS' , 9x, ' Area' ) 
12, 6x,Fl0.6, 6X,F10 .6, 6X,F10.5) 


TOTWE*16 


Write (6,737) 
strtot=0 
strcount=0 
f ortot=0 
f orcount=0 


Do 822 J*l, nmembs 
KA=start ( J) 

KB-finish (J) 

if (STRESS (KA, KB) .ge. 0) then 

Strrat-STRESS (KA, KB) /400 

strtot=strtot+Strrat 

strcount=strcount+l 

else 

Strrat=l 

endif 

If (FORCE (KA, KB) .It. 0) then 

Forrat=FORCE (KA, KB) /Per ( J) 
fortot-fortot+Forrat 
f orcount=f or count +1 

else 

Forrat=l 

endif 

WRITE (6,747) J, Forrat , Strrat , Area ( J) 

822 Continue 

Write (6^ *) 'Average safety factor in compression is' , -fortot /f orcount 
Write (6, *) 'Average safety factor in tension is' , strtot/streount 

Write (6, *)' That s all folks ' 

STOP 


ooo ooo ooo oooo 


END 


C****************Subroutine **************** ************ 
SUBROUTINE KMAT (C, THETA, KSPRING, I , J, nnodes) 


THIS TAKES K AS INPUT AND ADDS CONTRIBUTION TO 
CURRENT MEMBER 


REAL C ( 125, 125 ) , R (2 , 125 ) , KE (2 , 2 ) ,T(2,125) , KSPRING 


BEGIN BY FINDING R 

Na=2*nnodes 

DO 210 Ik-1,2 
DO 210 -11-1, Na 
210 R (Ik, II) -0 

R(l, 2*J-1) —COS (THETA) 

R(1,2*J)=S IN (THETA) 

R(2, 2*1-1) —COS (THETA) 

R (2 , 2*1) -SIN (THETA) 

DEFINE KE 

KE (1,1) -KSPRING 
KE (2, 2) -KSPRING 
KE(1, 2)— KSPRING 
KE (2, 1) —KSPRING 

MULTIPLY KE * R 

DO 220 11=1,2 
DO 220 JJ=1, Na 
T (II, JJ) —0 . 

DO 220 KK-1 , 2 

T (II, JJ) — T ( II , JJ) +KE ( II , KK) *R(KK, JJ) 
220 Continue 

DO 230 11=1, Na 
DO 230 JJ=l,Na 
DO 230 KK=1, 2 

C (II, JJ) =C ( II , JJ) +R (KK, II) *T (KK, JJ) 
230 Continue 

RETURN 

END 



1-1 


Appendix I 
Economic Analysis 


UPC= 


OCPF* 


MCPF- 


FCPF= 


NFLEET= 


NCYC= 


DESVOU 


DESWGTa 


DRANGE- 


FLVPD= 


FFPD= 


B 


376000 


600 


800 


321 


88001 


294001 


8 4 1 


ECONOMIC RESULTS 


UVCPF 


UWCPFb 


FLC ($ 


FCPOF ($fOZ 


FCPCI (S/CU. IN 


FLIFE {DAYS 


FCPVOUS/CU. IN. 


0.00032169 


0.00804214 


57069600 


20160000 

806400 


70.7708333 


2.83083333 


300 


6.47047619 


IN ORDER FOR PROPER PRICING SCHEMES. 


AVERAGE VALUES WERE USED NOT 


MAXIMUM VALUES 


CARGO LOAD PER DAY PER CITY 




M 


N 


TOTAL 


2060 


2760 


29101 


2620 


1560 


2450 


3450 


3950 


2830 


2310 


2500 1 


29400 


2450 


34501 


4050 
__27 30 
2210 


2600 


29 


EEUI 


M 

7001 

1 700 

N 

750 

l Z£Q 

TOTAL 

12420 

| 1 2420 


900| 

900 

155 

200 

i 


208 


1980 

161 0| 

1510 


1 6980 


iNON-OVERSEAS PRICE (S/CU, IN.) = | 8.7404762 


OVERSEAS PRICE (SFLIFE) = I 11.0104761 





















































































